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AN  ALTERNATIVE  APPROACH 
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State  space  methods  have  served  an  important  role  in  the  develop- 
ment of  modern  control  theory.  They  were  originally  proposed  to  help  solve 
the  increasingly  complex  aerospace  control  problems  which  frequency  response 
methods  of  the  time  could  not  deal  with.  Historically  this  theory  has  lead 
to  a much  greater  understanding  of  the  general  nature  of  a system  and  to  the 
solution  of  a great  variety  of  control  problems.  These  means  have  further 
been  applied  to  many  other  areas  in  the  sciences,  thus  firmly  establishing 
the  importance  of  modern  system  studies. 

While  the  elegance  of  these  methods  cannot  be  denied,  for  the 
linear  and  time  invariant  system  they  pose  many  practical  difficulties. 

These  include  incorporation  of  simplicity,  low  sensitivity,  and  robustness 
of  the  resulting  control,  as  well  as  integrity  and  fault  tolerance,  in  a 
design.  Working  with  an  input -output  description  of  the  plant  is  believed 
to  solve  many  of  these  problems  [1]. 

In  fact,  the  classical  frequency  domain  design  of  single  input 
single  output  systems  has  never  been  replaced  by  state  space  methods  in 
practice . 

Rosenbrock  et  al.,  at  the  Control  Systems  Centre,  University  of 
Manchester,  England,  has  extended  these  frequency  domain  techniques  to 
multiple  input  multiple  output  systems.  This  design  method,  which  he  has 
named  the  inverse  Nyquist  array  method,  encompasses  a large  number  of  control 
problems.  One  of  its  fortes  is  the  ability  of  the  designer  to  shape  the 

f 

outcome  of  the  resulting  designs  to  suit  his  needs  and  still  achieve  an 
adequate  system  performance. 


2 

How  this  is  possible  will  become  evident  in  what  follows. 

Chapter  2 relies  heavily  on  Rosenbrock's  own  exposition  of  the  method  [2] . 
Chapter  3 presents  an  outline  of  programs  to  be  used  for  interactive 
frequency  domain  design  with  computer  graphics  facilities.  An  example  is 
given  in  Chapter  4. 
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CHAPTER  2 

THE  INVERSE  NYQUIST  ARRAY  METHOD 


3 


2.1.  Introducing  Specifications  and  the  Plant 

The  design  process  begins  with  an  input-output  description  of  the 

plant.  Let  the  plant  be  linear  and  time  invariant.  It  has  an  (ixl)  vector 

u as  input,  and  an  (mXl)  vector  y as  output.  In  the  design,  one  is  primarily 

interested  in  the  response  at  each  output  of  the  plant  as  we  apply  sinusoidal 

forcing  functions  of  various  frequencies  to  each  input.  Therefore  the  plant 

is  represented  by  an  mxi  transfer  function  matrix  G(s),  which  is  rational, 

with  the  exception  of  the  inclusion  of  pure  exponential  terms  of  the  form 
sT 

e . These  terms  represent  pure  time  delay  and  may  be  included  when 
necessary,  as  many  plants  exhibit  this  sort  of  behavior. 

It  is  possible  that  a state  space  representation  of  the  plant, 
arising  from  a physical  model  or  available  empirical  data,  is  at  hand.  This, 
of  course,  is  a suitable  description  of  the  plant  after  the  required 
transfer  functions  are  obtained  via  analysis.  As  one's  actual  interest  is 
in  the  complex  frequency  mapping  between  input  and  output,  it  may  be  more 
appropriate  to  measure  this  data  on  the  act  ’>il  plant,  if  the  physical 
situation  allows  this.  For  the  type  plants  being  considered  the  data  may 
always  be  fitted  into  such  a transfer  function  matrix  expression.  While  this 
is  not  implicitly  a necessary  task  to  perform,  it  shall  be  assumed  that  the 
transfer  function  matrix  is  of  the  above  form.  In  the  case  of  designing 
from  actual  measured  data,  G(s)  merely  represents  this  data  and  is  a 
convenient  means  to  interpolate  between  the  data  points  for  the  physical 
model  at  hand. 
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As  usual  in  frequency  donialn  techniques,  one  deals  with  the  map  of 
the  open  loop  transfer  function  in  the  complex  plane  as  a suitably  defined 
contour  is  traversed.  The  encirclements  of  a critical  point  then  indicate 
stability.  Here  the  contour,  which  shall  always  be  denoted  by  D,  is  in 
every  case  traversed  clockwise.  It  consists  of  that  segment  of  the  imaginary 
axis  from  -jR  to  +jR,  as  well  as  the  semicircle  of  radius  R about  the  origin 
in  the  right  half  plane.  The  contour  is  indented  to  the  left  of  any  pole 
or  zero  of  the  transfer  function  at  the  origin  or  on  the  jtu  axis,  as  one 
in  general  is  interested  in  asymptotically  stable  systems  and,  therefore, 
must  include  these  points  within  the  closed  path.  The  radius  R is  taken  to 
be  large  enough  to  insure  that  all  right  half  plane  poles  and  zeros  are 
enclosed  by  D.  See  Figure  2.1.1. 

Given  a G(s),  consider  a design  procedure  based  on  the  closed 
loop  configuration  shown  in  Figure  2,1.2.  An  (£xk)  input  compensation  matrix 
K(s)  relates  the  error  E with  input  U.  Similarly  an  (kXm)  output  compensator 
matrix  L(s)  relates  output  Y to  z,  while  a G^Xk)  feedback  matrix  F(s)  is 
inserted  in  the  return  path. 

For  notational  simplicity,  the  (kxk)  open  loop  transfer  function 
matrix  of  the  system  relating  the  error  E(s)  to  the  output  Z(s)  is  defined 
as  Q 

G^  K?!).  (2.1.1) 

As  one  is  interested  in  the  closed  loop  responses,  the  (kXk)  closed  loop 
transfer  function  matrix  relating  inputs  V(s)  to  outputs  Z(s)  is  defined  as 
H(s).  By  assuming 


I 


Ik  + Q(s)  F(s)l  i 0 


it  immediately  follows  that 


H(s)  = (Ik+Q(s)  F(s))"  Q(s), 


(2.1.3) 


In  assuming 


L + F(s)  Q(s)|  ^ 0 


(2.1.4) 


the  closed  loop  transfer  function  may  be  expressed  as 


H(s)  Q(s)(Tj^+F(s)  Q(s))' 


(2.1.5) 


As  both  (2.1.3)  and  (2.1,5)  show  that  H is  a function  of  Q,  F,  and 


the  Laplace  variable  s,  one  might  consider 


H = hTq.F.s] 


(2.1.6) 


as  a compact  way  of  expressing  this  closed  loop  transfer  function  matrix. 

While  this  represents  a relatively  simple  expression  in  terms  of  computational 
considerations,  it  shows  very  little  in  terms  of  relating  the  open  loop 
system  to  the  closed  loop  system.  As  our  interest  is  in  designing  for 
an  acceptable  closed  loop  system  performance,  a further  investigation  into 
a possible  means  to  these  ends. becomes  necessary. 


2.2.  Inverse  Relationships  and  their  Consequences 

The  results  in  the  previous  section  are  valid  provided  that  the 


determinants  (2.1.2)  and  (2.1.4)  exist.  It  should  now  be  assumed  that 

IqI  * 0 (2.2.1) 
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which  implies 

1k|  ^ 0,  |g(  ^ 0 |l|  ^ 0.  (2,2.2) 

If  this  were  not  the  case,  the  system  would  not  be  functionally  controllable 
and  therefore,  be  of  very  little  practical  interest.  Functional  control- 
lability, also  known  as  functional  reproducibility,  loosely  means  that 
there  exists  some  input  which  will  generate  a specific  bounded  output  from 
zero  initial  conditions,  for  all  time  other  than  the  initial  time  [3,  p,  169]. 

Provided  that  (2.2.1)  is  satisfied,  Q will  have  an  inverse  which 
shall  be  denoted  by  a circumflex. 

= Q = (LGK)'^  = KGL  (2.2.3) 


A similar  notation  shall  be  used  for  all  other  matrix  inverses.  Note  that 
the  (i,j)  element  of  Q is  referred  to  as  whereas  the  inverse  of  that 
same  element  in  Q is  referred  to  as 

Using  inverse  matrices,  the  closed  loop  transfer  function  matrix 
implicitly  defined  in  (2.1.6)  has  an  Inverse  of 


H = F +Q. 


(2.2.4) 


This  shows  that  if  Q exists,  the  closed  loop  inverse  transfer 
function  can  be  expressed  in  the  astoundingly  simple  form  of  (2,2.4). 

Being  additive  by  nature,  it  is  the  convenient  means  of  relating  the  closed 
loop  system  to  the  open  loop  system  which  is  desired.  Clearly,  more 

insight  can  be  found  in  working  with  the  Inverse  system  and  the  easy 
transition  in  closing  the  loops,  than  from  working  with  the  noninversed 
system  and  dealing  with  the  nonintuitive  nature  of  the  matrix  multipliers 
imbedded  in  (2.1.6). 
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In  a design  one  is  interested  in  controlling  the  k outputs  Z with 
the  respective  k inputs  V.  To  do  this,  the  open  loop  inverse  transfer 
function  Q is  constructed  so  as  to  minimize  the  effect  of  the  off  diagonal 
entries  in  G.  Once  a sufficient  degree  of  noninteraction  is  present,  the 
original  problem  has  been  considerably  reduced. 

It  is,  of  course,  always  possible  to  construct  Q so  that  the 
desired  noninteraction  is  present  simply  by  choosing  k or  L equal  to  G. 

This  will  make  Q c agonal  and  equal  to  the  identity  matrix.  While  such  a 
method  has  been  suggested,  it  Imposes  severe  problems  in  terms  of  plant 
variations,  implementation  in  terms  of  stable  subsystems,  and  the  resulting 
shape  of  the  frequency  loci  [4] . Perhaps  the  most  objectionable  reason 
is  the  necessary  complexity  of  the  resulting  control  scheme.  As  simplicity 
in  a design  often  makes  the  difference  between  an  engineering  success  and  a 
failure,  one  could  never  justify  such  a scheme  as  this  with  its  inherent 
faults . 

Quite  the  contrary,  one  in  fact  wishes  to  find  a matrix  compensator 
which  is  as  simple  as  possible  in  order  to  achieve  a required  amcint  of 
noninteractiveness.  Ideally,  this  compensating  matrix  would  be  independent 
of  frequency  and  as  sparse  as  possible.  This  is  the  easiest  control  to 
implement  on  the  actual  plant,  and  would  probably  be  the  most  reliable  as 
well 

When  sufficient  decoupling  of  the  open  loop  system  is  attained. 


appropriate  stability  regions  can  be  determined  so  that  feedback  may  be 
applied.  By  considering 
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F = diagCf^]  , i«L,2,...,k  (2.2,5) 

where  each  is  scalar  and  independent  of  s,  the  closed  loop  frequency  loci 
of  the  diagonal  transfer  functions  in  H become  those  of  the  corresponding 
loci  in  f)  with  the  imaginary  axis  appropriately  shifted.  Single  loop 
frequency  compensation  in  the  classical  sense  of  phase-lag  and  phase-lead 
networks  may  be  designed  to  modify  these  frequency  loci  for  the  desired 
input-output  responses.  In  short,  our  multivariable  control  problem  becomes 
k single  loop  designs,  in  which  an  abundance  of  frequency  domain  techniques 
exist  [2,  pp.  28-116] , and  more  modern  algorithmic  methods  have  never 
entirely  replaced. 

Before  showing  how  this  is  possible,  it  is  necessary  to  formalize 
the  intuitive  idea  of  noninteractiveness.  As  the  design  means  adopted  are 
graphical  in  nature,  this  information  must  be  included  on  the  diagonal 
inverse  frequency  loci  of  each  individual  control  loop. 


2.3.  Dominant  Systems 

Rosenbrock  calls  these  loosely  coupled  or  noninteractive  systems 
dominant,  as  the  principle  diagonal  element  primarily  determines  that  input- 
output  relationship. 

2.3.1.  Diagonal  Dominance 

Definition  1 [2,  p.  142] : A complex  rational  (kxk)  matrix  Z(s)  is  said  to  be 

diagonally  row  dominant  on  the  contour  D if  Zj^j^(s)  has  no  pole  on  D for  all 


r 
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i * 1,2, . . , ,k,  and 


Zi^(s)|-  Silzij(s)|  > 0 


Vi*l,2,...,k 

and  for  all  s on  D. 


(2.3.1) 


Definition  2 [2,  p.  142]  : A complex  rational  (kXk)  matrix  Z(s)  is  said  to 


be  diagonally  column  dominant  on  the  contour  D if  Zj^^(s)  has  no  pole  on  D for 


all  i = 1,2, . . . ,k,  and 


lzii(s)|  -j|i|zji(s)|  >0  Vi»l,2,...,k 

and  all  s on  D. 


(2.3.2) 


Definition  3 [2,  p.  142]  : Z(s)  is  diagonally  dominant  on  D if  it  is  either 

diagonally  row  dominant  or  diagonally  column  dominant. 

Observe  that  row  dominance  implies  dominance,  and  column  dominance 
implies  dominance.  Then  for  dominance,  the  sum  of  the  moduli  of  the  off 
diagonal  row  elements  or  column  elements  must  be  greater  than  the  modules  of 
that  diagonal  element  for  all  s on  D.  It  is  clear  that  rows  and  columns 
may  not  be  mixed  to  check  dominance  at  a particular  s on  D.  Instead,  one 
may  have  either  row  dominance,  or  column  dominance  as  s takes  on  different 
values  around  D,  and  still  satisfy  Definition  3. 

A simple  graphical  construction  exists  to  check  the  dominance  of 
Q(s).  Let  q^^(s)  map  the  familiar  Nyquist  contour  D into  F^.  For  a 
particular  s on  D,  draw  a circle  with  its  center  at  the  corresponding  point 


on  of  radius 


di(s)  - j^iUij(s)l 


(2.3.3) 
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which  is  based  on  rows,  or 


d[(s) 


jjqji(s)| 


(2.3.4) 


which  is  based  on  columns.  Do  this  for  a sufficient  number  of  s on  D so 
as  to  define  an  area  about  F^.  If  the  area  defined  by  this  set  of  circles 
of  radius  d^  excludes  the  origin,  for  all  i*l,2,...,k,  then  Q is  row 
dominant  on  D.  Similarly  if  the  area  defined  by  the  set  of  circles  of 
radius  dl  excludes  the  origin  for  all  i-l,2,...,k,  then  Q is  column  dominant 
on  D.  In  this  same  sense,  if  none  of  the  circles  of  either  radius  d^  or  of 
radius  d^  has  the  origin  as  an  interior  point  or  on  its  circumference,  for 
all  i as  before,  then  Q is  dominant  on  D.  These  bands,  as  defined  by  (2.3.2) 
and  (2.3.3)  are  called  Gershgorin  bands  as  their  usefulness  is  dependent 
upon  Gershgorin  s Theorem. 


2.3.2.  Ostrowski's  Theorem 

Systems  which  are  diagonally  dominant  have  determinants  which 
show  interesting  properties  and  will  be  of  use  later  in  investigating 
stability.  However,  the  following  theorem,  due  to  Ostrowski  [5],  which  was 
rediscovered  by  Rosenbrock,  allows  for  the  location  of  the  diagonal  transfer 
function  in  a dominant  system. 


Gershgorin 's  Theorem,  while  being  of  little  practical  value  for 
this  design  method,  has  Important  consequences  in  the  rigorous  development 
of  what  follows.  For  its  statement  and  proof,  one  is  referred  to  the 
literature  [3,  p.  11] . 


I 
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Theorem  1:  Ostrowski's  Theorem  [2,p.  149]:  Let  the  rational  complex  (kXk) 

matrix  z(s)  be  row  (or  respectively  column)  dominant  for  s ■ on  any  simple 
closed  contour  C,  having  on  it  no  pole  of  any  z^^(s).  Then  Z(s^)  has  an 
inverse  which,  for  all  i“l,2,...,k  satisfies: 

^ (2.3.5) 

where 

0i(So>  * max  dj(SQ)/|zjj(sQ)l  (2.3.6) 

and  (^^^(Sq)  is  given  by  (2.3.3)  for  row  dominance.,  or  alternatively,  which 
satisfies 

^ ■‘I('o)  (2.3.7) 

where 

0^(Sq)  “ max  d](sQ)/lzjj(sQ)l  (2.3.8) 

jJi 

and  <1|(Sq)  is  given  by  (2.3.4),  for  column  dominance. 

The  proof  of  Ostrowski's  theorem  may  be  found  elsewhere  [3,pp.  204- 
206] . This  theorem  of  course  has  several  important  design  implications. 

Theorem  2 [2,  p.  150] : Let  0 exist  and  satisfy  the  conditions  of  Ostrowski's 
Theorem.  Furthermore,  assume  F is  diagonal  and  Independent  of  s.  Then 
for  each  s on  D the  diagonal  element  h^^^  of  (2.1.6)  satisfies 

lh'J(s)-(f^+q^^(s))l<  0^(s)d^(s)  < d^(s)  (2.3.9) 

if  H is  column  dominant,  or 

lh"J(s)-(fi+qii(s))l  < 0^(s)d^(s)  < d’(s)  (2.3.10) 


if  H is  row  dominant. 
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The  proof  of  this  theorem  follows  as  a direct  consequence  of 
Ostrovski's  Theorem.  Consider  Z ■ H * F+Q,  which  makes  ■ h^^  and 

^i  " ^i  “ ^i-^^ii- 


d.(s  ) 
1 ^ o^ 


0.  (s  ) = max  ^ 


(2.3.11) 


if  H is  column  dominant,  or 


0.'  (s  ) ■ max  i— 2 


(2.3.12) 


if  H is  row  dominant,  and  the  desired  result  is  obtained.  The  set  of 
circles  defined  by  (2.3.11)  and  (2.3.12)  lie  inside  the  corresponding 
Gershgorin  bands.  They  are  called  Ostrowski  bands  in  honor  of  the 
mathematician  on  whose  theorem  they  are  based.  Theorem  2 has  the  following 
graphical  consequences.  Pick  a particular  input-output  relationship 
characterized  by  a diagonal  element  in  H,  say  fi,,.  If  the  gains 


{ f , f2  > . . • , f j _ » f j ^2^  f • • • » 


(2.3.13) 


are  kept  constant  while  allowing  f^  to  vary,  then  one  has  essentially  a 
single  input  single  output  situation.  It  is  exactly  this  cause  and  effect 
relationship  for  which  we  wish  to  design  a single  loop  controller.  Then 
Theorem  2 tells  us  that  the  inverse  transfer  function  hj^^  lies  somewhere 
within  the  set  of  Ostrowski  bands  based  on  the  constant  gains  (2.3.13). 
Furthermore,  one  may  vary  fj  and  thus  design  for  the  appropriate  stability 
in  terms  of  gain  margins,  phase  margins,  other  Nyquist  diagram  type 
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specifications  [2,  pp.  33-35]  and  their  corresponding  time  domain  interpre- 
tations, from  the  inverse  Nyquist  diagram  with  the  appropriate  Ostrovski 
bands.  When  these  bands  are  small  enough,  the  resulting  composite  diagram 
is  treated  as  if  it  were  just  a "fuzzy"  inverse  Nyquist  diagram 


2.4.  Obtaining  Dominance 

In  Figure  2.1.2  the  forward  path  transfer  function  matrix  Q was 
forced  to  be  square  to  allow  the  possible  existence  of  its  necessary 
inverse.  This  was  done  by  making  matrices  K and  L of  the  appropriate 
dimensions.  It  now  becomes  necessary  to  make  the  further  assumption  that 
the  plant  transfer  function  matrix  G is  square.  This  implies  that  K and  L 
are  also  square  or  k = i*m. 

As  an  aside,  it  should  be  mentioned  that  one  does  not  wish  to 
limit  the  scope  of  this  method  by  assuming  a plant  with  an  equal  number  of 
inputs  and  outputs.  Rather,  at  the  present  time  there  is  no  general  way  of 
obtaining  a square  Q from  a nonsquare  G.  Clearly  more  work  is  needed  in 
this  area. 

When  confronted  with  a nonsquare  plant,  one  does  have  the  option 
of  reverting  to  the  direct  Nyquist  array  method  in  which  analogs  results 
as  to  dominance,  stability  and  the  location  of  the  closed  loop  transfer 
function  exist  [2,  pp.  174-179.  One  usually  chooses  not  to  pursue  this 
recourse  for  several  Important  reasons  [2,  p.  155]. 

A 

1)  Practical  experience  seems  to  show  that  there  is  a tendency  for  G 


to  be  more  dominant  than  G. 


I 

1 

I 

I 


2)  When  working  with  inverse  systems,  the  simplicity  of  relating  the 
open  loop  system  properties  to  those  of  the  closed  loop  system  adds 
much  additional  insight  into  a design. 

3)  Feedback  causes  the  width  of  the  Ostrowski  bands  to  rapidly  become 

small  in  the  inverse  system.  That  is  to  say,  as  feedback  is 

'.-I 

increased  in  the  other  loops,  we  make  h^  tend  towards  the 

uninversed  system  has  the  opposite  property  and  requires  all  the 
other  loop  gains  to  go  to  zero. 

As  our  primary  interest  is  in  designing  for  adequate  closed  loop  responses, 
it  is  clear  why  this  limitation  is  accepted  and  why  one  is  restricted  to  the 
inverse  system  representation, 

2.4.1.  Matrix  Operations 

Consider  the  (mXm)  input  compensator  matrix  K(s).  For  practical 
reasons  K(s)  should  have  all  its  poles  and  zeros  in  the  open  left  half 
plane.  This  is  becai  se  it  must  be  physically  implemented  by  a minimum  phase, 
asymptotically  stable  subsystem.  It  can  be  shown  [3,  p.  209]  that  such  a 
matrix  may  be  constructed  from 

= K^l^(s)k.^(s).  (2.4.1) 

When  working  with  inverses , this  is 

ft(s)  - &^(s)k^^(s)K^  (2.4.2) 

where  each  of  the  subscripted  square  matrices  has  the  following  function: 

1)  is  a permutation  matrix,  which  is  to  say,  in  each  column  and  each 
row  all  entries  are  zero  except  one,  which  is  unity.  Its  effect  is 


I 
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to  transpose  the  rows  of  any  matrix  which  it  premultiplies.  When 
operating  on  G its  function  is  to  renumber  the  inputs, 

2)  Kjj(s)  is  constructed  as  the  product  of  elementary  matrices  "^(s), 
each  of  which  has  unity  determinant  and  is  of  the  following  form: 


1 0 


0 1 


0 0 


-zf>(s)  . 


(2.4.3) 


0 0 


where  each  (s)  may  occur  in  any  off  diagonal  position  and  should 

be  a rational  function  of  s,  with  strictly  left  half  plane  poles 

^ (X) 

and  zeros.  The  effect  of  each  “^(s)  is  to  subtract  from  row  i 
a multiple  z^^  (s)  of  row  j from  the  matrix  which  premultiplies  it. 
When  operating  on  the  matrix  K^G(s),  the  effect  of  the  matrix  ^(s) 
is  to  recombine  the  parallel  paths  through  the  system  such  that  the 
resulting  matrix 


^^(s)k^G(s)  = ^n^(I^^^^(s))k^G(s) 


(2.4.4) 


is  diagonally  dominant. 

3)  The  matrix  K^(s)  is  diagonal  and  nonsingular.  As  it  represents  m 
single  loop  controllers,  which  will  be  designed  from  classical 
frequency  loci  techniques  once  sufficient  dominance  is  obtained 


from  (2,4.4),  K^(s)  should  contain  strictly  left  half  plane  poles 
and  zeros , 
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In  quite  an  analogous  fashion,  one  may  construct  the  (irt<m)  output 
compensator  matrix  L(s)  as 

L(iy  = ) La  (2.4.5) 

so  that 

L(s)  = LaLj^(s)L^(s)  (2.4.6) 

where  the  subscripts  indicate  matrices  whose  operations  correspond  to  those 
in  the  previous  case. 

Usually  one  designs  an  input  compensator  with  K(s)  as  in  (2.4.2) 
and  L = It  is  possible  to  design  an  output  compensator  with  K * and 

A 

L as  in  (2.4.6).  This  allows  some  additional  freedom  in  design  as  the 
physical  situation  requires.  Very  often  the  matrix  L(s)  will  be  moved  around 
the  loop  and  actually  implemented  in  the  input. 

One  may  always  attempt  to  construct  a diagonally  dominant  Q from 
these  elementary  row  or  column  operations.  Unfortunately,  this  becomes 
quite  a chore  for  large  systems.  Even  with  a plant  of  small  dimension  it 
is  very  possible  that  a satisfactory  solution  will  be  overlooked.  What  is 
needed  is  a more  systematic  method  of  obtaining  dominance. 

The  first  such  method  which  has  been  tried  is  to  diagonalize  G(o) 
with  a constant  k[6].  If  G(o)  is  nonsingular,  this  procedure  will  yield 
a Q(s)  which  is  dominant  about  w=o  , but  one  would  be  extremely  lucky  if 
dominance  is  maintained  throughout  the  D contour.  Similarly,  another 

A 

possibility  is  to  diagonalize  G(s)  as  s takes  on  large  imaginary  values, 
but  this,  too,  would  just  be  a guess.  Most  likely  dominance  would  fail 


elsewhere  on  the  contour. 
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One  could  attempt  to  combine  these  two  efforts  into  a more 
complicated  expression,  such  as 

K(s)  = D^K^+k^s  (^.7) 

where  diagonalizes  G(o),  diagonalizes  G(s)  as  s approaches  j»,  and 
is  a diagonal  matrix.  Unfortunately  for  systems  of  many  dimensions,  the 
necessary  inversion  of  K(s)  given  by  (2.4.7)  becomes  rather  complicated, 
is  awkward  to  calculate  exactly,  and  would  be  difficult  to  implement  on 
the  actual  physical  plant.  Some  of  these  problems  may  be  overcome  if, 
instead,  one  used  the  form 

K(s)  = D k +D,K, /s  (2.4.8) 

'•^ooll 

where  diagonalizes  G(s')  as  s approaches  j «> , diagonalizes  g7o),  and 
and  are  each  diagonal  (gain)  matrices. 

While  either  of  the.  forms  (2.4.7)  or  (2.4.8)  may  be  tried,  the 
difficulty  is  obtaining  a dominant  Q at  intermediate  frequencies  along  the 
Wyquist  D contour  remains.  Furthermore,  our  interest  is  still  in  finding 
a compensating  matrix  which  is  as  simple  as  possible.  A satisfactory 
constant  compensating  matrix  might  exist  which  would  suite  the  purpose. 

A 

This  leads  one  to  investigate  possible  diagonalization  of  G(s)  at 
jo  < s < j»  , Of  course,  the  dominance  of  the  resulting  Q will  have  to  be 
checked  on  the  entire  domain  of  s in  question,  but  a method  which  gives 
such  a possible  candidate  is  needed. 


2.4.2.  Pseudodiagonf.llzatlon  [7] 

The  diagonalization  of  Q(juJ)  by  choosing 
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k = G(jcu) 


(2.4.9) 


would  require  K to  be  a function  of  s,  as  G is  generally  a complex  function. 
One  could  assume  K to  be  constant  and  formulate  the  problem  as  minimizing 
the  sum  of  the  moduli  of  the  off  diagonal  terms  of  the  resulting  0(s).  Upon 
placing  some  convenient  constraint  on  K,  it  can  be  shown  that  the  resulting 
minimization  problem  is  straightforward  to  formulate,  but  impossible  to 
solve  for  a strictly  real  matrix  K [ 7] . 

This  problem  may  be  overcome  if  one  tries  instead  to  minimize 
the  sum  of  the  squares  of  the  moduli  of  the  off  diagonal  terms.  While  this 
may  be  done  at  any  specific  s = jtu,  instead  the  problem  is  generated  at  N 
discrete,  independently  weighted  intervals  of  frequency  to^..  That  is,  by 
considering  row  j of  Q,  the  problem  is  to  find  row  j of  k such  that  the 
function 


N 

S 

r=l 


r 

^kSl-'V 

k5^j 


(2.4.10) 


is  minimized.  Here  each  Y^.  is  chosen  to  be  positive  and  real,  and  represents 
the  weighting  factor  for  each  ou^  as  r=l,2,...,N.  While  the  following 
results  are  based  on  an  input  compensator  and  rows  as  this  is  the  usual 
case,  analogous  results  for  output  compensators  or  columns  may  be  obtained 
in  a similar  manner. 


Element 


9-i,(j“^  ) expression  (2.4.10)  may  be  written  as 
jK  r 


m 


(2.4.11) 
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and  while  is  complex,  it  is  expressed  by  real  and  imaginary  parts  as 

2lkO»r>  -“i'k'+Jelk'- 


Then  the  function  (2.4.10)  to  minimize  becomes 


liy  E|.Efc..(c(5^jp|-)|23. 

r=l  rk=l'i=l  rk  ■^'^ik  ' 

k?^j 


(2.4.13) 


Upon  introducing  the  convenient  constraint  of 


m 


Efc:,  = 1 

i=l  Ji 


(2.4.14) 


and  appending  the  constraint  with  the  Lagrange  multiplier  X,  expression 
(2.4.10)  requires  the  minimization  of  the  function 


m , m 


k?tj 


(2.4.15) 


Taking  the  necessary  partial  derivatives  with  respect  to  k^^  and  setting 
them  equal  to  zero  to  obtain  the  required  minimum  results  in 


I 

I 

I 

I 


Sk 


kJ^j 


For  simplicity  the  row  vector  Kj  is  defined  as 


itj  ■ (fcjp 


Similarly  the  matrix 


^ r-l^rl-k-l^^ik  ®ik  ^^k  ^£k 
kJ^j 


(2.4.16) 


(2.4.17) 


(2.4.18) 


is  introduced,  Aj  is  symmetric  and  real,  and  at  least  positive  semidefinite 
by  construction,  so  all  of  its  eigenvalues  are  real  and  nonnegative.  This 
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allows  (2.4.16)  to  be  written  in  the  form  of  a standard  eigenvalue  problem 


A.feT  -XkT  = 0. 
J J J 


(2.4.19) 


While  any  eigenvector  of  Aj  satisfies  expression  (2.4.19),  if  one  constructs 
kj  from  the  smallest  nonzero  eigenvalue  of  A^,  the  minimum  of  the  original 
function  (2.4.10)  is  obtained. 

This  procedure  yields  a jth  row  of  k which  is  of  unit  length  and 
minimizes  the  sum  of  the  squares  of  magnitude  of  the  off  diagonal  terms  of 
row  j in  Q(jou).  It  is  clear  that  this  must  be  done  for  all  the  rows 
j = l,2,...,m.  Nothing  prevents  the  use  of  different  values  of  or 
from  (2.4.10),  or  just  a single  frequency  and  unity  weighting  for  all  the 

A 

rows.  In  any  case,  the  diagonal  dominance  of  the  resulting  Q(s)  must  always 
be  checked  at  all  points  of  interest. 

The  constraint  (2.4.14)  leads  to  a relatively  easy  computational 
problem.  Unfortunately,  while  satisfying  (2.4.14)  it  often  will  simul- 
taneously lead  to  unacceptably  small  values  of  To  overcome  this. 


the  alternate  constraint 


Al5jjO”r>l  ■ ‘ 


(2.4.20) 


may  be  substituted  for  (2.4.14).  A similar  analysis  shows  that  the  more 
difficult  solution  of  the  eigenvector  problem 


(2.4.21) 


is  required.  Here  Aj  is  as  in  (2.4.18),  and  Ej  is  the  symmetric,  positive 
semidefinite  matrix  defined  by 
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^ r=l^r'-"ij  +Pij  J- 


(2.4.22) 


The  idea  of  pseudodiagonalization  may  be  extended  by  considering 
compensating  matrices  dependent  upon  frequency  as  in  (2.4.7)  where  higher 
power  terms  in  s may  be  included  as  needed.  The  general  comments  of  the 
previous  subsection  still  motivate  one  against  such  a choice.  The  problem 
is  further  compounded  in  that  minimization  techniques  as  this  in  no  way 
guarantees  a resulting  matrix  K(s)  which  has  strictly  left  half  plane  poles 
and  zeros . 

As  can  be  seen,  obtaining  open  loop  dominance  is  very  often  a 
function  of  the  physical  situation  at  hand.  The  degree  of  dominance  is 
usually  dependent  upon  the  extent  of  complexity  the  system  designer  will 
allow  to  be  actually  implemented.  Certainly  before  these  methods  can  be 
put  into  a more  general  use,  additional  work  on  obtaining  dominance  is 
needed. 


2.4.3.  Closed  Loop  Dominance  via  Dynamic  F(s) 

One  may  consider  obtaining  dominance  in  the  closed  loop  system. 

If  G(o)  is  nonsingular,  there  always  exists  a constant  input  compensating 
matrix  K as 

K - G(o)  (2.4.23) 

because  G(o)  is  real.  This  will  make  Q(o)  diagonal  and  dominance  will 
remain  for  at  least  a small  region  close  to  the  origin  in  the  complex  plane. 
Consider 


* 
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F(s)  = F(o)  +F^(s)  (2.4.24) 

where  F(o)  is  diagonal  and  scalar.  Fj^(s)  has  zero  diagonal  elements.  Its 
off  diagonal  elements,  which  are  a function  of  frequency,  are  chosen  so  that 
the  resulting  dominance  of  H in  the  closed  loop  system  is  improved.  As 
Fj^(o)  is  constructed  to  be  zero,  the  resulting  F(s)  yields  transient 
compensation  of  the  interaction. 

The  scope  of  the  method  has  not  yet  been  well  defined  [8] . Such 
a method  has  never  actually  been  physically  implemented.  If  it  would  be, 
great  care  in  the  amount  of  phase  advance  introduced  and  its  effect  on  the 
dynamic  behavior  of  the  system  should  be  carefully  considered.  The 
sensitivity  of  the  system  with  respect  to  changes  in  the  plant  will  also 
have  to  be  investigated,  as  a small  change  in  Q may  severely  change  the 
degree  of  dominance  when  cancellation  of  terms  in  F+Q  are  large. 


2.5.  Stability 

In  multivariable  control  problems,  the  allowable  gain  in  one  loop 
is  very  often  a complicated  function  of  the  gains  in  the  other  loops.  For 
design  purposes  one  would  like  to  have  a good  idea  of  what  combinations  of 
gains  result  in  asymptotic  closed  loop  behavior.  Permissable  feedback  values 
could  then  be  chosen,  perhaps  with  the  additional  requisite  that  stability 
is  maintained  as  these  gains  vary  a prescribed  amount  from  their  design 
values.  In  addition,  it  may  be  desirable  to  design  the  closed  loop  system 
so  that  it  does  not  become  unstable  in  the  case  of  some  catastrophic  event. 


I 


such  as  a transducer  or  actuator  failure.  Before  showing  how  this  is 
possible  some  results  of  a general  nature  are  necessary. 


2.5.1,  Frequency  Domain  Criteria  for  Stability 

Theorem  3 [2,  p.  141];  Let  Q have  an  inverse  Q as  defined  by  (2.2.3),  which 

has  p^  poles  in  the  right  half  plane.  Let  |q|  map  the  Nyquist  contour  D into 

f and  IhI  map  D into  f . As  s goes  once  around  D,  will  encircle  the 

y H Q 

A ^ '' 

origin  N times,  while  F encircles  the  origin  times,  where  all  positive 
Q H H 

encirclements  are  clockwise.  Then  the  resulting  closed  loop  system  is 
asymptotically  stable  if  and  only  if 


Nq-\=  Po- 


(2.5.1) 


The  proof  of  this  theorem  requires  showing  that  the  return 
difference  matrix  relates  the  open  loop  poles  to  the  closed  loop  poles  and 
must  be  motivated  on  a very  general  ground  from  the  theory  of  polynomial 
system  matrices  and  the  "Principle  of  the  Argument"  in  complex  number  theory. 
As  it  is  not  the  exact  result  we  are  seeking  for  design  purposes,  its  proof 
may  be  found  elsewhere  [2,  pp.  131-141], 

T^  eorem  4 [2,  p,  143] (Dominant  Encirclement  Theorem):  Let  the  (kXk)  rational 

complex  matrix  Z(s)  be  dominant  on  any  simple  closed  contour  C which  has 
on  it  no  pole  of  z^^(s)  for  all  i = l,2,...,k.  Let  z^^j^(s)  map  C into  F^,  which 
encircles  the  origin  times,  and  lz(s)j  map  C into  F^  , which  encircles  the 
origin  times.  All  clockwise  encirclements  are  again  taken  to  be  positive. 


N = , 

z i=l  i 


(2.5.2) 
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The  proof  of  the  Dominant  Encirclement  Theorem  relies  on  Gershgorln's 
Theorem  and  a necessary  corollary.  It  is  rather  lengthy  and  may  be  found 
elsewhere  [2,  pp.  25-26]. 


From  these  two  theorem  statements,  a number  of  different  stability 
theorems  may  be  motivated  depending  upon  which  matrices  are  assumed  to  be 
dominant,  and  whether  or  not  one  is  working  with  inverse  systems.  Rather 
than  stating  all  of  the  possibilities,  one  chooses  instead  to  give  a result 
which  is  the  most  useful  one  for  design  purposes. 

Theorem  5:  (Generalized  Graphical  Inverse  Nyquist  Stability  Criterion)* 

Let  Q and  H be  dominant  on  D,  and  assume  that  F is  diagonal  and  independent 


of  s . 


Suppose  each  of  the  Gershgorin  bands  based  on  the  diagonal 
elements  of  Q exclude  the  origin  and  the  critical  point  (-f^,0). 
Let  these  bands  encircle  the  origin  N . times  and  encircle  the  critical 
point  (-fj^,0)  times.  Then  the  resulting  closed  loop  system  is 

asymptotically  stable  if  and  only  if 


^ A 

.r  N . - S N 
i=l  qt  i=l  hi 


(2.5.3) 


where  p is  the  same  as  in  Theorem  3. 

Proof : When  F is  diagonal  and  independent  of  s,  exclusion  of  the  origin  by 

A A 

a set  of  bands  based  upon  h^^^  = ^i^^ii  same  as  the  exclusion  of  the 

point  (-f^,0)  by  a set  of  bands  based  on  While  noting  that  the  premise 

of  the  Dominant  Encirclement  Theorem  applies  and  gives  the  total  number  of 


Additional  theorems  and  a more  thorough  development  leading  up  to 
the  statement  of  this  theorem  may  be  found  in  Rosenbrock  [2,  pp.  143-154]. 
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required  encirclements  in  both  the  open  and  closed  loop  system,  Theorem  3 is 
used  and  the  required  results  are  obtained. 

This  theorem  represents  the  Generalized  Inverse  Nyquist  Stability 
Criterion  for  multivariable  systems  and  has  the  graphical  design  implications 
as  stated.  It  should  be  mentioned  that  it  is  possible  to  exclude  the  origin 
with  one  set  of  bands,  for  example  based  on  rows,  while  excluding  the 
critical  point  (-f^,0)  by  a set  of  bands  based  on  columns,  as  is  otherwise 
implicitly  obvious  in  the  use  of  the  word  dominance.  While  this  theorem  is 
of  most  use  in  its  present  form,  it  sometimes  is  possible  to  make  H dominant 
but  not  Q.  In  this  special  case,  there  is  the  following  result 
Corollary  1 [2,  p.  145] : Let  H be  dominant  on  D and  satisfy  the  dominant 

encirclement  theorem.  Then  the  closed  loop  system  is  asymptotically  stable 
if  and  only  if 


■ Pq  - "-Q  - »o 


(2.5.4) 


where  S(,  . and  p are  as  in  Theorem  5 and  p„  and  are  the  poles  and  zeros 
hi  o Q Q 

of  1q1. 

Sketch  of  Proof:  (2.5.1)  in  Theorem  3 states  that  for  asymptotic  stability, 

one  must  have 


N - N = P 
Q H ‘^o 


(2.5.5) 


while  the  dominant  encirclement  theorem  states 


\ ■ i=l\i’ 


(2.5.6) 


An  application  of  the  Principle  of  the  Argument  gives 


V ‘’q  ■ 


(2.5.7) 
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and  combining  (2.5.5)  through  (2.5.7)  yields  the  desired  statement.  Notice 
also  that  these  results  can  be  useful  if  |q|  is  known  algebraically. 

Underlying  the  stability  theorems  in  either  the  open  or  closed 
loop  case  is  the  idea  of  diagonal  dominance.  If  a system  cannot  be  made 
dominant,  the  theorems  say  nothing. 


2.6.  A Finishing  Remark 


If  any  design  method  is  to  be  of  good  use  it  must  be  physically 
practical,  both  in  terms  of  implementation  and  measurement.  In  particular. 


the  inverse  Nyquist  array  method  fails  in  this  respect  as  those  values  of 
complex  frequency  which  correspond  to  the  semicircular  arc  on  the  Nyquist 
contour  D have  no  physical  significance. 

There  is  a natural  desire  to  be  relieved  from  considering  this 
portion  of  the  Nyquist  contour  in  a design.  Not  only  is  it  impossible  to 
perform  actual  measurements  on  the  system  for  these  values  of  s,  dominance 
may  fail  in  this  region  as  well. 

While  theorems  have  been  presented  which  lead  one  to  believe  this 
is  possible,  after  Rosenbrock,  the  following  is  stated  as  an  assumption 
[2,  pp.  153-155] : 

'Assumption:  Those  poles  which  determine  stability  cross  the 

segment  of  the  imaginary  axis  between  s = -JuJq  and  s = where 

oUq  is  some  specified  (real)  frequency.' 

Therefore,  one  can  investigate  stability  by  considering  imaginary  values  of 

s from  0 < s < Here  is  greater  than  and  > R,  as  defined  by 


the  inclusion  of  all  right  half  plane  poles  and  zeros  of  the  transfer 
functions  by  the  Nyquist  contour  D. 

Being  able  to  exclude  this  semicircular  arc  makes  the  practical 
considerations  of  the  design  method  complete. 


I 
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CHAPTER  3 

COMPUTER  PROGRAMMING 

3.1.  Introduction 

Rosenbrock's  method,  as  outlined  in  the  previous  chapter, 
represents  a feasible  alternative  in  the  design  of  multivariable  control 
systems.  It  does  require  a rather  large  amount  of  calculation,  even  for 
plants  of  low  order.  This  could  be  done  by  hand  and  the  results  graphed 
accordingly,  but  the  necessary  time  and  labor  would  make  such  an  approach 
extremely  unappealing. 

Fortunately,  with  the  present  availability  of  high  speed,  time 
sharing  computers  this  computational  effort  can  be  completely  eliminated. 
Furthermore,  with  the  advent  of  graphical  display  terminals,  these  results 
may  be  immediately  graphed  and  examined.  Finally,  with  the  development  of 
interactive  software,  a suitable  control  scheme  can  be  built  up  as  insight 
gained  into  the  system  is  quickly  reapplied  until  an  acceptable  design  is 
reached , 

Such  an  effort  has  been  undertaken  using  a Tektronix  graphics 
terminal  with  supporting  software  on  the  Coordinated  Science  Laboratory's 
DEC- 10  computer. 

Programming  in  the  form  of  twenty-one  subroutines  has  been 
developed  in  DEC  system- 10  FORTRAN  code  [9]  and  may  be  compiled  using  either 
the  FlO  or  the  F40  FORTRAN  compiler.  Graphic  routines  have  been  written 
with  the  aid  of  the  Tektronix  Advanced  Graphics  AG-II  supportive  software 
package  [ 10] . 


F 


mmmmmm 


31 

All  the  subroutines  are  functionally  oriented  and  have  been  kept 
general,  to  be  useful  in  the  design  of  any  control  problem.  The  routines 
operate  on  data  which  is  stored  in  a common  block. 

As  conversational  software  for  the  analysis  and  design  of  systems 
has  become  in  vogue  [11],  it  is  hoped  that  these  subroutines  will  be 
included  as  an  option  package  in  an  interactive  software  program.  At  that 
time  the  common  data  base  will  have  to  be  appropriately  modified  to  conform 
with  the  conventions  already  established  in  the  particular  program  used. 

It  is  possible  to  combine  the  written  subroutines  themselves  into 
an  interactive  software  program.  An  example  of  how  this  may  be  done  is 
included  with  the  next  chapter  in  investigating  a particular  design 
example . 

3.2.  Common  Block  Usage 

The  common  data  area  has  been  divided  into  nine  separately  labeled 
segments,  each  of  which  serves  its  own  functional  importance.  For  efficient 
computer  utilization  when  solving  a specific  problem,  these  common  areas 
should  be  dimensioned  as  the  particular  situation  dictates.  The  routines 
are  therefore  stored  in  separate  files  without  the  common  data  base  and 
commented  as  to  which  labeled  common  block  area  is  required. 

COMMON  /Al/  Block 

Labeled  common  area  /Al/  stores  all  of  the  system  matrices  as 
shown  in  the  block  diagram  structure  included  in  Figure  2.1.2.  Matrices 
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K(s)  and  L(s)  are  each  broken  into  constant  matrices,  XKA  and  XLA,  and  into 
frequency  dependent  matrices,  XKB  and  XLB,  respectively. 

The  transfer  function  of  the  plant  is  stored  in  the  four 
dimensional  array  G as  follows:  considering  location 

G(I,J,K  L)  (3.2.1) 

defines  a storage  position  where 
I = row  index  of  transfer  function 
J = column  index  of  transfer  function 
K = coefficient  of  (K-1)  power  of  s in  G(I,J) 

L = 1 for  the  numerator  polynomial 
= 2 for  the  denominator  polynomial. 

Then  when  dimensioning  G,  it  is  clear  that  I and  J must  be  equal 
and  correspond  to  the  number  of  inputs  or  outputs,  K must  be  one  more  than 
the  highest  power  of  s occurring  in  any  entry  of  G,  and  L is  2. 

The  coefficients  of  the  polynomial  terms  in  matrices  XKB  and  XLB 
are  referred  to  in  exactly  the  same  manner  and  are  dimensioned  as  needed. 

Matrices  XKA,  XLA,  and  F are  all  square,  constant,  and  of 
appropriate  dimension. 

COMMON  /A2/  Block 

Labeled  common  area  /A2/  stores  seven  single  variables.  In  order 
of  occurrence  they  are: 

OMEG0  - real  variable  denoting  the  starting  frequency  of  the  Nyquist  path. 


I 

I 

I 


It  will  almost  always  be  defined  to  be  zero,  but  is  included  so  that 
a region  on  the  ju)  axis  not  including  the  origin  may  be  investigated. 


33 


OMEGND  - real  variable  corresponding  to  the  end  frequency  on  the  Nyquist 
path  on  the  ju)  axis.  Referring  to  Figure  2.1.1,  it  is  also  the 
length  of  the  radius  R defining  the  semicircular  arc  in  the  right 
half  plane. 

DTOMEG  - real  variable  corresponding  to  the  change  in  frequency  from  OMEG0 
to  OMEGND  as  one  travels  up  the  jou  axis. 

ISTART  - integer  variable  denoting  the  starting  location  in  arrays  yet  to 

be  defined.  While  it  usually  is  1,  it  has  been  included  so  we  may 
plot  any  group  of  points  corresponding  to  a particular  segment  of 
the  frequency  loci, 

lEND  - integer  variable;  it  defines  an  end  location  as  ISTART  defines  a 
starting  location.  On  the  jou  axis,  one  is  considering 


TEND  = [ (OMENGD-OMEG0) /DTOMEG]  + 1 


(3.2.2) 


total  number  of  such  points. 

NP  - integer  variable  denoting  the  number  of  points  between  each 
Gershgorin  or  Ostrowski  band  plotted. 

MO  - integer  variable  correspond  to  number  of  inputs  or  outputs.  In 
reference  to  Chapter  2,  it  is  the  order  of  the  system  or  k, 

COMMON  /A3/  Block 

Labeled  common  area  /A3/  is  a scratch  area  used  as  temporary 
storage  in  many  of  the  computational  routines.  It  has  been  further  divided 
into  three  separate  areas. 

Area  /A31/.  Stores  two  square  complex  matrices  Cl  and  C2,  each  which  should 


be  dimensioned  (MOX  MO).  It  also  has  two  real  variable  vectors  Si  and  S2  of 


length  MO. 
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Area  /A32/.  Stores  two  real  variables,  three  dimensional  arrays  XRl  and 
XII,  which  correspond  to  real  and  imaginary  parts  of  the  Nyquist  map.  Each 
should  be  dimensional  as  (MOx  MOx  lEND).  When  referring  to  the  point  in 
the  complex  plane  stored  in  position 

(XR1(I,J,N),XI1(I,J,N))  (3.2.3) 

then  one  is  considering  the  (I,J)  input-output  relationship  at  a frequency 
corresponding  to  the  Nth  point  on  the  frequency  loci.  For  points  on  the  jou 
axis , the  frequency  is 

CD  = (N-1)*DT0MEG  (3.2.4) 

as  N^will  take  on  values  from  ISTART  to  lEND. 

COMMON  A/33/.  Stores  two  real  variables,  three  dimensional  arrays  XR2  and 
XI2,  each  dimensioned  as  (MO x MO X lEND) . They  define  positions  and  serve 
a function  similar  to  arrays  XRl  and  XII  in  area  /A32/. 

COMMON  /A4/  Block 

Labeled  common  area  /a4/  stores  the  real  and  imaginary  parts  of 
the  open  loop  inverse  transfer  function  Q.  It  does  so  in  two,  three 
dimensional  real  variable  arrays,  QR  and  QI  respectively.  Each  is  dimen- 
sioned as  (MOx  MOX  lEND)  and  the  (I,J,N)  storage  position  contains  informa- 
tion in  the  same  manner  as  defined  around  (3.2.3)  and  (3.2.4). 

COMMON  /A5/  Block 

Labeled  common  area  /A5/  stores  information  including  the  value  of 
the  Gershgorin  bands  in  the  four  dimensional  real  value  array  BAND.  Stored 


in  location 


BAND  (I,J,N,L) 


is  the  sum  of  the  magnitudes  not  including  the  (I,J)  magnitude  at  frequency 
N of  the  Q values  from  area  /A4/,  The  sum  of  the  row  magnitudes  are  stored 
in  L = 1 and  the  sum  of  the  column  magnitudes  are  stored  in  L*2.  Array  BAND 
should  be  dimensioned  as  BAND  (MO, MO, TEND, 2) . 

COMMON  /A6/  Block 

Labeled  common  area  /A6/  contains  one  three  dimensional  real 
valued  array  OST,  which  stores  the  values  of  the  Ostrowski  bands.  In 
location  OST  (I,N,L)  is  the  Ostrowski  band  based  on  h^^^  at  frequency  N where 
again  L=1  for  the  band  based  on  rows  and  L = 2 for  the  band  based  on  columns. 
Array  OST  should  therefore  be  dimensioned  as  OST  (MO, TEND, 2). 

COMMON  /a7/  Block 

Labeled  common  area  /A7/  contains  two  important  real  value  arrays 
XS  and  XC,  XS  is  square  and  dimensioned  as  XS(MO,MO).  XC  is  four  dimen- 
sional and  dimensioned  the  same  as  G in  area  /Al/  as  defined  by  (3.2.1).  In 
the  subroutines  which  follow,  calculations  are  performed  with  matrices  XS 
and  XC  which  will  take,  on  values  from  area  /Al/  and  thus  allow  the  design 
of  either  input  or  output  compensators. 


3.3.  Subroutines 

The  subroutines  have  been  divided  into  three  separate  sections 
which  in  some  respect  correspond  to  the  operation  which  they  perform.  Each 
individual  routine  is  stored  separately  on  disk  under  a file  name  which 
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corresponds  to  the  subroutine  name  with  an  extension  of  SAV.  Common  block 
usage  other  than  COM^K)N  /A2/  has  been  appropriately  commented  but  not 
included.  Subroutines  have  further  been  commented  as  to  what  interactive 
procedures  they  perform,  if  any,  and  what  additional  subroutines  are  called 
in  that  routine,  if  any,  A listing  of  the  subroutines  and  a brief  description 
of  each  of  their  functions  is  included  in  Appendix  2. 

3.3.1.  Library  COMP 

Library  COMP  contains  seven  computational  subroutines  useful  in 
performing  the  necessary  calculations  at  various  stages  of  a design. 

Subroutine  MAP 

The  subroutine  MAP  computes  the  complex  value  of  the  transfer 
function  matrix  G(s)  as  s takes  on  values  around  the  Nyquist  path  as  defined 
by  the  real  variables  in  COMMON  /A2/.  \ 

The  calling  sequence  is 

CALL  MAP  (ICIRCL) 

where 

ICIRCL  = 0 maps  the  value  of  G(s)  for  strictly  imaginary  values  of  s on 
the  ju)  axis  from  OMEG0  to  OMEGND  in  discrete  steps  of  DTOMEG. 

The  total  number  of  these  points  is 

NTIMES  = ((OMEGND-OMEG0 ) /DTOMEG)  + 1.  (3.3.1) 

ICIRCL  = 1 maps  the  value  of  G(s)  as  s goes  around  the  semicircle  of  radius 


OMEGND,  It  does  this  in  99  discrete  steps. 
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ICIRCL  = 2 (or  any  other  integer)  causes  s to  accept  values  along  the 
entire  Nyquist  path.  It  will  do  so  in 


(NTIMES  +99  +NTIMES) 
discrete  steps. 


(3.3.2) 


The  resulting  value  of  G(s)  is  a complex  number  whose  real  part 
is  stored  in  XRl  and  whose  imaginary  part  is  stored  in  XII  in  COMMON  /A32/, 
for  each  g^^^  and  the  N such  discrete  steps  performed. 

If  a pole  of  G is  encountered  anywhwere  on  the  imaginary  axis, 
subroutine  MAP  types  on  the  screen  the  message 


A POLE  AT  OMEGA  = XXXXX.XX  FOR  G(X,X,XXX). 

It  then  adjusts  the  value  of  the  complex  frequency  variable  s accordingly  to 

s = (OMEGA,  -DTOMEG)  (3,3.3) 

so  that  these  poles  are  included  inside  the  contour  D. 

Subroutine  MAG 

The  subroutine  MAG  takes  the  complex  value  of  the  functions 
implicitly  defined  in  array  XC  and  stores  them  in  the  labeled  common  area 
/A32/  in  arrays  XRl  and  XII.  It  does  so  at  discrete  frequencies  defined  by 
(3.2.4)  where 

N * 1,2 lEND.  (3,3,4) 


This  subroutine  is  used  in  the  design  of  frequency  dependent 
compensators.  The  appropriate  array  from  COMMON  /Al/  is  loaded  into  XC  and 
MAG  computes  ita  complex  value.  Note  that  if  one  were  to  set  XC = G,  then 
MAG  performs  the  same  function  as  calling  MAP(0). 
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Singularities  of  functions  in  XC  have  been  programmed  around  in 
the  same  fashion  as  in  routine  MA.P.  The  user  is  notified  of  these  as  the 
statement 

A POLE  AT  OMEGA  = XXXXX.XX  FOR  FUNCTION  (X,X,XXX) 

is  typed  on  the  terminal  screen. 

The  calling  sequence  is 

CALL  MAG 

and  has  no  calling  arguments. 

Subroutine  MIMP 

The  subroutine  MIMP  performs  the  matrix  multiplication  of  constant 
and  matrix  XS  from  COMMON  /A7/  with  the  complex  matrix  represented  by 
(XRljXIl)  from  COMMON  /A32/.  It  does  so  for  the  N positions  defined  by 
(3,3,4).  The  resulting  complex  matrix  is  appropriately  stored  in  labeled 
common  /A4/  area  by  its  real  part  in  array  QR  and  its  imaginary  part  in 
array  QI. 

Routine  MIMP  finds  its  usefulness  in  the  design  of  constant  matrix 
compensators.  It  also  presents  an  easy  means  of  transferring  data  from 
labeled  common  area  /A32/  to  labeled  common  area  /A4/  when  matrix  XS  is 
set  equal  to  the  identity  matrix  I^.  As  all  the  graphics  routines  plot 
frequency  loci  from  data  contained  in  arrays  QR  and  QI  from  common  /A4/, 
this  is  a necessary  transfer  to  perform  in  persuing  the  map  of  G(s). 

The  calling  sequence  is 

CALL  MIMP 

where  no  calling  arguments  are  needed. 
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Subroutine  MCMP 

The  subroutine  MCMP  is  similar  to  routine  MIMP  except  that  it 
performs  the  multiplication  of  two  complex  matrices.  The  routine  takes  the 
product  of  the  complex  matrix  defined  by  XRl  and  XII  from  COMMON  /A32/  with 
the  complex  matrix  defined  by  XR2  and  XI2  from  COMMON  /A33/.  It  then  stores 
the  resulting  matrix,  again  by  real  and  imaginary  parts  in  QR  and  QI  in 
COMMON  /A4/.  This  is  performed  for  all  N as  in  (3.3.4). 

The  routine  is  useful  in  the  design  of  frequency  dependent 
compensators.  The  calling  sequence,  which  has  no  arguments,  is 

CALL  MCMP 

Subroutine  GERSH 

The  subroutine  GERSH  computes  the  moduli  of  the  Q values  from 
COMNK)N  /A4/  and  stores  the  resulting  real  number  in  array  BAND  in  COMMON 
/A5/  as 

MO 

BAND(I,J,N,1)  = (QR(I,K,N),qi(I,K,N))l  (3.3.5) 

MO  I , 

BAND(I,J,N,2)  = Z 1 (QR(K,J,N),QI(K,J,N))1 . (3.3.6) 

K*  X 

kjtl 

When  I = J the  value  stored  in  BAND(I, J,N,1)  represents  the  radius 
d^  from  expression  (2.2.3)  for  Gershgorin  bands  based  on  the  row  s of  Q, 
while  the  value  stored  in  BAND(1, J,N,2)  represents  the  radius  d^  from 

A 

expression  (2.3.4)  for  Gershgorin  bands  based  on  the  columns  of  Q.  This  is 
for  all  N as  in  (3.3.4).  Note  that  those  values  in  BAND  for  I J are  not 
actually  Gershgorin  bands  but  the  appropriate  sum  of  the  magnitudes  including 
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the  diagonal  term.  This  information  may  be  useful  in  determining  how  a 
preliminary  input  compensator  might  be  constructed  in  order  that  the 
diagonal  dominance  is  achieved,  and  therefore  has  been  included. 


Upon  executing  the  statement 


CALL  GERSH 


these  magnitudes  are  computed  from  QR  and  QI  and  stored  in  BAND. 


Subroutine  OSTROW 

The  subroutine  OSTROW  computes  the  Ostrovski  bands  and  places 
the  resulting  value  of  the  circular  radius  in  array  OST  as  discussed  in 
the  labeled  common  area  /A6/  section.  These  computed  values,  as  defined 
by  (2.3.11)  and  (2.3.12),  are  based  on  the  ^ values  from  COMION  /A4/, 
their  resulting  diagonal  BAND  values  from  COMMON  /A5/  and  the  current  feed- 
back gains  in  diagonal  matrix  XS  from  COMMON  /a7/.  In  the  closed  loop 
inversed  system,  as  one  changes  the  amount  of  feedback,  the  size  of  the 
Ostrovski  bands  change.  It  then  becomes  necessary  only  to  call  subroutine 
OSTROW  for  each  change  of  feedback,  as  the  Q values  and  BAND  values  remain 
constant  as  the  loops  are  closed. 

The  Ostrovski  bands  are  calculated  by  executing  the  statement 


CALL  OSTROW. 


Subrout ine  MMIJ 


The  subroutine  MMIJ  determines  the  minimum  and  maximum  points  on 
the  composite  plot  of  a Nyquist  diagram  with  its  Gershgorin  bands.  It 
therefore  requires  use  of  the  Q values  in  COMMON  /A4/  and  the  BAND  values 
in  COMMON  /A5/. 
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The  calling  sequence  is 

CALL  J,XMIN,XMAX,YMIN,YMAX) 

where 

I = row  index  of  composite  plot 
J = column  index  of  composite  plot 
XMIN  = minimum  horizontal  value  of  composite  plot 
XMAX  = maximum  horizontal  value  of  composite  plot 
YMIN  = minimum  vertical  value  of  composite  plot 
YMAX  = maximum  vertical  value  of  composite  plot. 

Subroutine  MMIJ  searches  the  composite  diagram  for  I = J and 
compares  the  MIN  and  MAX  variables  with  the  outer  points  of  the  composite 
diagram.  For  I ^ J the  same  search  is  performed  on  the  frequency  loci  without 
the  Gershgorin  bands.  In  either  case  it  is  necessary  to  initialize  the  MIN 
variables  to  a large  positive  value  and  the  MAX  values  to  a large  negative 
value  in  calling  the  routine. 

This  subroutine  serves  two  purposes  when  used  in  conjunction  with 
the  graphics  routines.  First,  it  determines  where  the  screen  margins  should 
be  set  to  when  drawing  a composite  diagram,  as  the  supportive  software  for 
the  graphics  terminal  has  no  way  of  doing  this  when  the  Gershgorin  bands  are 
included.  Second,  it  provides  a means  by  which  a common  scale  may  be 
determined  when  plotting  several  such  diagrams  on  the  screen  at  once. 

3.3.2.  Library  GRAPH 

Library  GRAPH  contains  six  subroutines  which  display  all  the 
information  necessary  in  performing  a design. 
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Subroutine  PTNYl 

The  subroutine  PTNYl  plots  the  Q values  in  COMMON  /A4/  from 
position  ISTART  to  position  TEND  on  the  terminal  screen.  The  supportive 
software  automatically  sets  screen  limits  by  the  CALL  CHECK  statement  if 
these  have  not  been  previously  set  elsewhere.  The  CALL  DSPLAY  statement 
causes  the  Nyquist  plot  and  axis  to  be  graphed  on  the  screen. 

The  calling  sequence  is 

CALL  PTNYl (I, J) 

where 

I = row  index  of  loci  plotted. 

J = column  index  of  loci  plotted. 

Subroutine  PTGRl 

The  subroutine  PTGRl  causes  the  Gershgorin  bands  to  be  graphed 
over  the  Nyquist  plot  by  use  of  the  supportive  software  routine  CPLOT. 

It  plots  these  bands  starting  from  position  ISTART,  and  plots  one  band 
every  NP  positions,  to  position  lEND. 

The  calling  sequence  is 

CALL  PTGRl (I, J, IRC) 

where 

I = row  index  of  bands  plotted 
J = column  index  of  bands  plotted 
IRC  = 1 for  Gershgorin  bands  based  on  rows 

* 2 for  Gershgorin  bands  based  on  columns. 


F 
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Subroutine  PTOSTl 

The  subroutine  PTOSTl  plots  the  Ostrovski  bands  in  exactly  an 
analogous  fashion  as  routine  PTGRl  plots  the  Gershgorin  bands. 

The  calling  sequence  is 

CALL  PTOSTl (I, IRC) 

where 

I ■ diagonal  position  of 
IRC  * 1 for  Ostrovski  bands  based  on  rows 

* 2 for  Ostrovski  bands  based  on  columns. 

Subroutine  PTBIGl 

The  subroutine  PTBIGl  interactively  plots  a single  Nyquist  plot 
and  the  corresponding  Ostrovski  bands  using  the  maximum  possible  terminal 
screen  size. 

Upon  executing  the  calling  statement 
1 CALL  PTBIGl 

I 

the  graphic  display  panel  is  erased  and  the  message 

^ TYPE  0 TO  RETURN,  OR  TYPE  1,1/2  FOR  H(1,I)  AND  OST  ROWS/COLUMNS 

is  typed  at  the  top  of  the  screen.  Entering  0 again  erases  the  screen  and 
returns  execution  to  the  calling  program.  If  values  are  entered,  the  screen 
is  cleared  and  the  supportive  software  is  adjusted  accordingly  so  that  the 
desired  composite  diagram  is  plotted  on  the  screen  using  routines  PTNYl 
and  PTOSTl.  The  same  message  as  before  is  then  typed  at  the  very  upper 
edge  of  the  screen,  allowing  the  process  to  repeat  itself  as  often  as 


needed . 
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Subroutine  MAK4 

The  subroutine  MAK4  causes  the  four  Nyquist  diagrams  corresponding 
to  the  upper  left  (2X2)  square  portion  of  Q to  be  graphed  on  the  screen. 

It  also  plots  the  Gershgorin  bands  for  the  two  diagonal  entries  of  the 
defined  submatrix  if  desired.  While  the  system  under  investigation  could 
be  constructed  so  that  this  routine  would  display  something  useful,  its 
most  meaningful  application  is  for  a two  input,  two  output  system. 

The  routine  first  searches  all  four  of  the  composite  diagrams  for 
the  maximum  and  minimum  horizontal  and  vertical  points,  by  four  successive 
calls  to  subroutine  MMIJ.  This  data  is  then  used  with  some  supportive 
software  routines  to  set  the  scale  of  each  plot  to  be  equal.  The  Nyquist 
diagrams  are  graphed  on  the  screen  with  subroutine  PTNYl  and  the  Gershgorin 
bands,  if  desired  are  plotted  on  the  diagonal  loci  with  routine  PT®1. 

The  calling  sequence  is 

CALL  MAK4  (IRC) 

where 

IRC  * 0 causes  no  Gershgorin  bands  to  be  plotted 

A 

* 1 causes  Gershgorin  bands  based  on  the  row  elements  in  Q to  be 
plotted 

A 

= 2 causes  the  Gershgorin  based  on  the  column  elements  in  Q to  be 
plotted . 

Subroutine  DRW 

Subroutine  DRW  is  the  major  interactive  graphics  routine.  It 
basically  serves  a routing  function  and  relies  on  the  previous  subroutines 
to  actually  perform  a given  operation.  As  such,  it  does  not  deal  with  any 


I 

» 


physical  data  and  does  not  require  any  common  block  usage.  Included 
within  subroutine  DRW  is  entry  point  DRWl.  A flow  chart  of  the  routine 
is  included  as  Figure  3.3.1. 
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Upon  executing  the  statement 

CALL  DRW 

the  message 

PLOT  DIAGRAM:  TYPE:  0 FOR  NO,  1 FOR  YES 

is  typed  on  the  terminal  screen.  Entering  0 causes  execution  to  return  to 
the  calling  program,  while  entering  1 causes  execution  to  proceed  to  entry 
point  DRWl.  This  entry  point  has  been  included  as  usually  one  determines 
whether  or  not  a graphical  output  is  desired  from  the  calling  program. 

If  this  is  the  case  and  it  has  been  previously  determined  that  a plot  is 
desired  one  may  instead  execute  the  statement 

CALL  DRWl  . 

In  either  case,  at  this  point  the  message 

ENTER  -0,0-  FOR  4 or  -I,J-  for  Q(I,J) 

is  typed  on  the  terminal  screen.  Entering  nonzero  integer  values  causes 
subroutine  MMIJ  to  search  for  the  necessary  MAX  and  MIN  values,  the  scaling 
is  set  and  subroutine  PTNYl  plots  the  appropriate  frequency  loci. 

If  two  zeros  are  entered,  or,  after  the  single  loci  is  graphed, 
the  message 

PLOT  BANDS?  TYPE  0 FOR  NO,  1 FOR  ROWS,  2 FOR  COLUMNS 


is  typed  on  the  terminal  screen.  In  the  case  of  plotting  four  diagrams, 
subroutine  MAK4  with  the  proper  band  identifying  argument  is  called  and 
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the  four  frequency  loci  with  the  desired  bands  are  graphed  on  the  terminal 
screen.  In  the  case  of  a single  plot,  if  Gershgorir  bands  are  to  be  plotted, 
subroutine  PTGRl  with  the  proper  band  identifier  is  called. 

In  either  case  the  routine  then  routes  itself  back  to  its  beginning 
and  the  original  message 

PLOT  DIAGRAMS?  TYPE  0 FOR  NO,  1 FOR  YES 

is  typed  on  the  screen  as  before,  allowing  the  process  to  repeat  itself 
as  many  times  as  desired. 


3.3.3.  Library  AUX 

Library  AUX  contains  eight  subroutines  of  an  auxiliary  nature 
which  perform  the  interactive  modification,  disk  storage  or  outputting  of 
various  variables. 

Subroutine  M0DIA2 

The  subroutine  M0DIA2  allows  the  three  integer  variables  ISTART, 
lEND,  and  NP  in  labeled  common  area  /A2/  to  be  modified  from  the  computer 
console  during  a design. 

i 

Upon  executing  the  calling  statement 


CALL  M0DIA2 


the  message 


ENTER  ISTART,  TEND,  NP,  OR  AND  RETURN 


is  typed  on  the  screen.  Entering  three  integer  values  appropriately 
modifies  these  variables,  while  a comma  for  any  entry  causes  them  to  remain 


I 

I 

I 

I 
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the  same.  The  routine  then  types  the  three  values  on  the  screen  in  a 
unformated  output  so  there  is  no  mistake  as  to  what  values  they  have 
become . 

Subroutine  MODXS 

The  subroutine  MODXS  performs  the  interactive  modification  or 
output  of  matrix  XS  in  COMMON  /A7/  as  diagramed  in  Figure  3.3.2. 

Upon  executing  the  calling  statement 


the  message 


CALL  MODXS 


MODIFY  XS?  TYPE  0 FOR  NO,  1 BOR  YES 


is  typed  on  the  screen.  Entering  1 causes  the  statement 
TYPE  I,J,X  IN  XS(I,J)  * X 

to  be  typed  on  the  screen.  After  the  values  are  entered,  the  first  message 

MODIFY  XS?  TYPE  0 FOR  NO,  1 FOR  YES 

again  appears,  allowing  additional  entries  in  XS  to  be  changed. 

Upon  returning  0 the  statement 


OUTPUT  XS?  TYPE:  0 FOR  NO,  1 FOR  TTY,  2 FOR  LP 


is  typed  on  the  screen.  After  entering  the  appropriate  value  the  output 
is  performed  and  execution  is  returned  to  the  calling  program. 

Subroutine  MODXC 

The  subroutine  MODXC  performs  the  interactive  modification  and 


output  of  array  XC  as  routine  MODXS  does  for  matrix  XS.  This  again  may  be 
represented  by  the  diagram  included  as  Figure  3.3.2. 


Enter 

Values 


Output 


TTY  Output 


LP  Output 


Return 


Functional  flow  diagram  for  subroutines 
MODXS  and  MODXC. 


Figure  3.3.2 
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The  calling  statement  is 

CALL  MODXC  . 

This  causes  the  statement 

MODIFY  XC?  TYPE  0 FOR  NO,  1 FOR  YES 

to  be  typed  on  the  terminal  screen.  When  returning  1 the  statement 

I,J,N,X,l/2  IN  XC(I,J)  = X*(S)**N,NUM/DOM 

is  typed  on  the  screen  and  after  the  values  are  entered,  the  original 
message  again  appears,  allowing  possible  further  modification  in  array  XC. 
Entering  0 causes  the  statement 

OUTPUT  XC?  TYPE  0 FOR  NO,  1 FOR  TTY,  2 FOR  LP 

to  appear  on  the  screen.  Returning  a value  causes  the  appropriate  output 
to  be  performed  and  execution  is  then  returned  to  the  calling  program. 
Subroutine  POUT 

The  subroutine  QOUT  causes  the  Q values  in  COMMON  /A4/  from 
positions  ISTART  to  lEND  to  be  dumped  for  numerical  investigation. 

The  calling  sequence  is 

CALL  QOUT(IDUM) 

where 

IDUM  * 1 for  terminal  screen  output 
IDUM  * 2 for  line  printer  output 
Subroutine  QDSKWR 

The  subroutine  QDSKWR  stores  all  the  ^ values  in  COMMON  /A4/  on 


disk  file  Ql.DAT. 
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Upon  executing  the  calling  statement 


CALL  QDSKWR 


the  message 

STORE  Q VALUES?  TYPE  0 FOR  NO,  1 FOR  YES 


is  typed  on  the  terminal  screen.  Returning  0 causes  a return  to  the  calling 
program.  Typing  1 causes  the  formated  storage  of  the  integer  1 followed  by 


I,J,N,QR(I,J,N),QI(I,J,N) 


on  the  disk  file  Ql.DAT.  These  values  may  be  perused  after  execution  by 
outputing  the  file,  or  used  for  the  next  program  run  by  using  the  subroutine 
QDSKRD . 

Subroutine  QDSKRD 

The  subroutine  QDSKRD  reads  the  Q values  from  disk.  It  assumes 
that  the  disk  file  Ql.DAT  exists  and  looks  at  the  value  stored  in  the  second 
space  of  its  first  line.  If  this  value  is  0,  a return  is  made  to  the  calling 
program  and  the  Q values  are  not  read.  If  instead  this  value  is  1,  the  Q 
values  are  read  and  stored  in  arrays  QR  and  Q1  in  COMMON  /A4/ . The  routine 
also  adjusts  the  variable  lEND  in  COMMON  /A2/  as  the  largest  Nth  position 
read  into  these  arrays. 

The  calling  sequence  is 


CALL  QDSKRD  (IDUM) 

where 

IDUM  ■ 0 on  return  indicating  Q values  were  not  read 

■ 1 on  return  indicating  Q values  were  read  and  stored. 
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Subroutine  OUT 

The  subroutine  OUT  outputs  all  the  system  matrices  stored  in 
C0M1K)N  /Al/  and  all  the  programming  variables  from  COMMON  /A2/  to  either  the 
terminal  screen  or  the  line  printer.  It  also  has  the  ability  to  call  sub- 
routine QOUT  for  a look  at  the  current  Q values  from  ISTART  to  lEND. 

Upon  executing  the  calling  statement 


CALL  OUT 


the  message 


FOR  OUTPUT  TYPE 

0 FOR  NONE,  1 FOR  TTY,  2 FOR  LP 


appears  on  the  terminal  screen.  Typing  0 causes  a return  to  the  calling 
program.  Typing  1 or  2 causes  the  message 


FOR  SYSTEM  OUTPUT  TYPE  0 
FOR  Q OUTPUT  TYPE  1. 

Entering  0 causes  the  system  and  programming  variables  to  be  dumped, 
whereas  typing  a 1 causes  the  subroutine  QOUT  to  be  called. 

In  either  case,  after  the  appropriate  output  is  performed,  the 

message 


TYPE  0 TO  CONTINUE,  TYPE  1,T0  DO  AGAIN 

to  appear  on  the  terminal  screen.  Entering  0 causes  execution  to  return  to 
the  calling  program,  while  entering  1 causes  the  original  message 

FOR  OUTPUT  TYPE 

0 FOR  NONE,  1 FOR  TTY,  2 FOR  LP 
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to  be  typed  on  the  screen.  This  allows  the  process  to  repeat  itself  again, 
if  desired. 

Subroutine  FINISH 


The  subroutine  FINISH  is  executed  by  the  calling  statement 


CALL  FINISH. 


The  message 


YOU  ARE  DONE : BEST  PRINTOUT  SYSTEM.'  I ! 


is  typed  on  the  terminal  screen  and  subroutine  OUT  is  called  so  that  the 
current  system  values  may  be  saved.  The  screen  is  then  erased  and  subroutine 
QDSKWR  is  called  so  that  Q values  may  be  stored  for  further  use.  Execution 
is  then  terminated. 


1 


1 
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4.  DESIGN  ORIENTATION  AND  EXAMPLE 


4.1.  Introduction 

The  subroutines  introduced  in  the  previous  chapter  have  been 
combined  into  a main  program  and  three  design  packages.  These  have  been 
applied  to  the  particular  design  problem  of  a lateral  autopilot  to  maintain 
leading  and  roll  attitude. 

The  programming  presented  here  is  oriented  to  the  two  input,  two 
output  autopilot  plant  using  input  compensators.  Output  compensators  may 
be  designed  in  an  analogous  fashion  by  modifying  the  'work'  matrices  from 
COMPTON  /A7/.  Particular  attention  should  be  given  to  the  order  of  matrix 
multiplication  for  these  types  of  designs.  Plants  of  higher  dimensions 
cause  no  problem  either,  if  the  common  data  area  is  appropriately  constructed 
as  was  discussed  in  Section  3.2. 


4.2.  Autopilot  Plant 

The  lateral  motions  of  an  aircraft  as  defined  by  the  roll,  yaw 
and  sideslip  angles  are  coupled  together,  but  are  almost  completely  uncoupled 
to  the  longitudinal  or  pitch  and  plunge  motions.  The  resulting  perturbation 
equations  to  first  order  of  the  fifth  order  system  describing  these  lateral 
motions  are  found  [12]  to  be: 
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(4.2.2) 

(4.2.3) 

(4.2.4) 

(4.2.5) 

State  variables  are  defined  as 

0 = roll  angle 
p = roll  angular  velocity 
t = yaw  angle 
r * yaw  angular  velocity 
0 = sideslip  angle 

and  inputs  are 

6r  = rudder  deflection 
6a  = aileron  deflection 

where  (4.2.2)  and  (4.2.3)  must  be  solved  simultaneously  for  a proper  state 
space  description.  Doing  so  and  using  typical  values  of  the  coefficients 
for  an  aircraft  weighing  100,000  pounds  flying  at  500  miles  per  hour  at 
30,000  feet  [12l,  the  resulting  state  space  representation  of  the  plant  is 
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(4.2.6) 


Computer  analysis  [ll]  shows  the  system  to  be  controllable  and 
have  eigenvalues  of 

« 0.0 

* --886 

(4.2.7) 

X^  = +.0045 

X,  , = +.026  ± j .644 

4,5 

indicating  an  open-loop-unstable  system. 

Before  a design  can  begin,  a suitable  choice  of  output  variables 
must  be  made.  Intuitively  one  would  control  the  rudder  deflection  with  the 
yaw  angle  and  the  aileron  deflection  with  the  roll  angle.  Investigating 
the  original  regulator  design  made  from  minimizing  a performance  index 
[12]  shows  this  primarily  to  be  the  case,  where  the  respective  angular 
velocities  also  show  high  gains  in  the  corresponding  feedback  path. 

In  light  of  this  information,  the  angles  *1'  and  0 are  chosen  as 
the  outputs  of  the  system.  These  angular  velocities  could  instead  have 


57 


been  chosen,  but  the  all  Integrator  block  diagram  of  the  plant  (Figure  4,2.1) 
shows  this  would  result  in  some  unobservable  states  and  would  regulate  the 
angular  velocities  and  not  the  angular  positions  of  the  aircraft. 

Computer  analysis  [ll]  shows  the  system  to  be  observable  and 
routinely  yields  a transfer  function  description  of  the  plant  as 


■[.3807  s^+.3121  s^-f.00893  S+.01951  -s[,0404  s^+.0503  s -f  .447] 


DETI 


L^2J 


[.06715  s^-.0178  s^ 000588  s -I-, 0226]  s[l,587  s^+. 0624  s + .602] 


6r 

6a 


(4.2.8) 


where  DET  “ s[s^  + .829  s^+.364  s^+. 365 s -.00166] 


(4.2.9) 


4.3.  Autopilot  Main  Program 

Before  the  main  program  for  the  autopilot  design  may  be  run,  two 
files  must  be  created  on  disk.  The  first  is  file  Ql.DAT,  which  is  used 
to  store  the  Q values  from  COMMON  /A4/,  if  desired.  The  integer  0 is  stored 
in  the  second  space  of  the  first  line,  and  the  file  is  left  with  no  line 
numbers.  This  procedure  is  illustrated  in  Appendix  1. 

The  second  file  which  must  be  created  is  the  data  file  DATA12.F0R. 
This  contains  the  data  subroutine  DATA12,  which  initializes  the  system  data 


Figure  4.2.1.  All  integrator  block  diagram  showing  state  space 
description  of  autopilot  plant  with  BLTPLTS. 
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in  COMMON  /Al/  and  the  programming  variables  in  COMMON  Ikll  via  standard 
FORTRAN  assignment  statements.  A copy  of  DATA12.F0R  is  included  in 
Appendix  2,  along  with  all  the  autopilot  programming  routines. 

The  main  program  is  executed  from  the  command  file  MAIN.CMD 
as  illustrated  in  Appendix  1.  Upon  its  execution,  the  graphics  terminal 
and  supportive  software  are  initialized,  lEND  is  calculated,  and  subroutine 
DATA12  is  called  for  the  necessary  data. 

The  subroutine  QDSKRD  is  then  called.  If  disk  file  Ql.DAT  has 
Q values  stored  in  it  from  the  previous  run,  these  values  are  automatically 
reloaded  into  arrays  QR  and  QI  in  COMMON  /A4/  and  the  routine  returns  the 
integer  value  IDUM  = 1.  This  causes  the  message 

Q VALUES  WERE  OBTAINED  FROM  DISK 

to  be  typed  on  the  terminal  screen  and  execution  proceeds  to  one  of  the 
six  main  program  options. 

If  no  Q values  are  stored  from  the.  previous  run,  QDSKRD  returns 
the  integer  value  IDUM  = 0.  This  causes  matrix  XS  to  be  set  equal  to  the 
Identity  matrix,  subroutine  MAP  (0)  is  called,  causing  the  value  of  G(s) 
along  the  jw  axis  to  be  stored  in  COMMON  /A32/.  Subroutine  MIMP  is  called, 
taking  these  values  and  placing  them  in  COMMON  /A4/  also. 

The  subroutine  QINV  is  then  called.  This  subroutine  takes  the 
inverse  of  the  (2x2)  square  complex  matrices,  lEND  in  number,  stored  in 
COMMON  /A4/.  It  returns  the  inversed  values  again  to  COMMON  /A4/.  Note 
that  one  could  have  started  with  the  inversed  system  data  in  routine  DATA12. 
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Taking  the  point  by  point  inverse  with  QINV  saves  one  the  labor  of  computing 
the  inverse  transfer  function  matrix  by  hand.  Subroutine  QINV  is  stored 
on  disk  file  QINV2.FOR  and  is  included  in  Appendix  2. 

In  either  case,  COMMON  /A4/  now  stores  the  open  loop  system 
frequency  data-  The  message 
TYPE:  0 TO  STOP 
TYPE:  1 FOR  OPEN  LOOP  PACKAGE 
TYPE:  2 FOR  INPUT  COMPENSATOR  PACKAGE 
TYPE:  3 FOR  CLOSED  LOOP  PACKAGE 
TYPE:  4 FOR  OUTPUT 

TYPE:  5 TO  MODIFY  INTEGER  COMMON  /A2/  VALUES 

is  typed  on  the  screen.  Entering  the  desired  Integer  causes  execution  to 
proceed  to  a particular  conversational  routine.  Returning  from  one  of  these 
routines  causes  this  same  message  to  appear  and  allows  the  user  to  chose 
a different  option,  until  he  decides  to  stop. 


4.4.  The  Interactive  Design  Packages 

Entering  0 to  stop  causes  subroutine  FINISH  to  be  called.  Similarly 
entering  4 calls  routine  OUT  and  entering  5 calls  routine  M0DIA2.  These 
subroutines  are  included  with  the  other  five  auxiliary  type  subroutines 
in  the  file  AUX  FOR. 

In  this  same  manner  the  six  computationed  routines  have  been 
copies  from  their  respective  .SAV  file  into  the  file  COMP. FOR.  The  seven 
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graphics  routines  are  similarly  stored  in  GRAPH  FOR.  The  necessary  common 
block  usage,  dimensioned  for  the  autopilot  example,  is  included  in  each 
subroutine  of  these  three  fi^es. 

4.4.1.  Subroutine  OLISP 

Typing  1 for  the  open  loop  package  calls  the  subroutine  OLISP. 

This  routine  is  useful  in  investigating  the  diagonal  dominance  of  the  plant 
and  in  the  design  of  constant  inptit  compensators.  It  can  also  be  used  to 
show  the  effect  of  actuator  failure  at  an  input  to  the  plant. 

The  routine  has  a functional  flow  as  diagrammed  in  Figure  4.4.1. 
Upon  entering  the  routine,  work  matrix  XS  is  set  equal  to  input  compensating 
matrix  XKA.  The  open  loop  system  Q variables  are  stored  in  arrays  XRl  and 
XII  in  COMMON  /A32/  and  subroutine  GERSH  computes  the  Gershgorin  bands  for 
those  values.  The  message 

DEALING  W/  OPEN  LOOP  INPUT  CONST.  COMP 


is  typed  on  the  screen  to  remind  the  user  of  the  action  being  taken.  The  j 

routing  message 

TYPE  0 TO  RETURN,  1 TO  PLOT,  2 TO  MODIFY 

i 


then  appears  on  the  terminal  screen.  Entering  2 calls  subroutine  MODXS  for 
modification  and/or  intermediate  output  of  the  constant  input  compensator 
and  upon  its  return  the  same  directioning  message  is  typed  on  the  screen. 
Entering  0 to  return  sets  system  matrix  XKA  equal  to  the  routine's  work 


matrix  XS,  and  causes  the  execution  to  continue  at  the  main  program. 


start 


XS  =XKA 
Save  Q Values 
in  Area/A32/ 


Compute  Gershgorin 
Bands  (GERSH) 


Plot  (1)  I , ' I Modify  (2) 

' Return/Plot/Modify 


XS  ^(XRl.XIl) 
Stored  in  Q 
(MIMP) 


Compute  Gershgorin 
Bands  (GERSH) 


Graphical  Output 
(DRW) 


1 

Modify /Output 
XKA  Values 
(MODXS) 

1 

1 

! 

f 

Set  XKA  = XS 


Return 


Figure  4.4.1.  Functional  flow  chart  of  subroutine  OLISP. 
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If  instead  one  chooses  to  display  a graphical  output  by  entering  1, 
the  following  happens.  Sub.^utine  MIMP  takes  the  matrix  product  of  matrix 
XS  and  the  original  Q,  values  stored  in  COMMON  /A32/  and  places  the  result 
in  arrays  QR  and  QI  in  COMMON  /a4/.  Subroutine  GERSH  computes  the  moduli 
of  the  new  Q values  and  subroutine  DRW  is  called  at  entry  point  DRWl  for  a j 

graphical  display  as  described  in  Section  3.3.2.  After  all  the  desired  plots 
are  Inspected,  the  routing  message 

TYPE:  0 TO  RETURN,  1 TO  PLOT,  2 TO  MODIFY 

again  appears  on  the  screen,  allowing  further  modification  and  plotting,  or  ' 

a return  to  the  main  program.  i 

When  XS  is  set  equal  to  the  identity  matrix,  the  information 
displayed  is  the  open  loop  plant.  Otherwise  it  is  the  compensated  plant. 

The  important  thing  to  note  is  that  when  returning  to  the  main 
program,  the  values  of  Q are  those  of  the  last  input  compensator  used.  If 
it  is  not  desired  to  change  these  values  from  what  they  were  when  the 
package  was  entered,  matrix  XS  should  be  set  equal  to  the  identity  matrix 
and  the  results  plotted  once.  The  effect  of  successive  changes  in  XS  are 
not  accumulative  unless  a return  is  made  to  the  main  program. 

One  may  exploit  this  particular  quirk  as  follows.  Let  the  plant 
Q values  be  stored  on  disk.  This  is  usually  the  case  at  this  stage  in  a 
design.  Enter  this  design  package  and  set  matrix  XS  equal  to  the  identity 
matrix,  except  for  one  diagonal  entry.  This  is  set  to  zero  to  represent 
an  actuator  failure  at  that  input.  Return  to  the  main  program,  then  return 
to  this  open  loop  design  package  and  set  XS  equal  to  the  desired  design 
value  if  the  input  compensator.  The  resulting  frequency  loci  will  be  those 


when  this  actuator  has  failed.  When  compensator  XKA  is  not  diagonal,  this 
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type  of  failure  may  have  astounding  effects  on  the  resulting  stability 
region  of  the  closed  loop  system.  Note  that  the  Q values  are  partially 
destroyed  by  this  type  of  procedure,  and  therefore  should  be  saved 
beforehand. 

4.4.2.  Subroutine  IFDCD 

Typing  2 for  the  input  compensator  package  option  in  the  main 
program  calls  subroutine  IFDCD.  This  routine  is  very  similar  to  the  open 
loop  package  OLISP.  The  difference  between  the  two  is  that  here  the  input 
compensator  is  allowed  to  be  a function  of  frequency.  It  therefore  finds 
its  usefulness  in  those  types  of  designs,  but  it  may  be  used  for  investi- 
gation of  the  plant  frequency  loci  or  show  the  effect  of  actuator  failures 
as  in  subroutine  OLISP.  A functional  flow  chart  of  routine  IFDCD  is 
included  as  Figure  4.4.2. 

Upon  entering  this  package,  the  systems  Q values  from  COMMON  /a4/ 
are  stored  in  arrays  XR2  and  XI2  in  labeled  common  area  /A33/.  The  work 
matrix  XC  is  set  equal  to  the  input  frequency  dependent  matrix  XKB  and 
subroutine  GERSH  is  called  to  compute  the  sum  of  the  moduli  of  the  Q values. 
The  message 

DEALING  W/  FREQ.  DEPNDT.  INPUT  COMP, 
is  typed  on  the  screen. 

The  user  then  has  the  same  three  options  as  in  the  previous 
routine  as  the  routing  message 


TYPE:  0 TO  RETURN,  1 YO  PLOT,  2 TO  MODIFY 


Save  Q Values 
in  Area/A33/ 


Compute  Gershgorin 
Bands  (GERSH) 


Plot  (1) 


Re  turn/ Plot /Modify 


Take  Complex 


Modify /Output 
XKB  Values 
(MODXC) 


Compute  Matrix 
Product  of 


XKB  * G 

Store  in  G (MCMP) 


Return  (0) 


Return 


Figure  4.4.2.  Functional  flow  chart  of  subroutine  TFDCP 
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Bands 
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> 

Graphical  Output 
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appears  on  the  terminal  screen.  Entering  2 calls  subroutine  MODXC  to  modify 
values  in  the  inverse  transfer  function  of  the  input  compensator  and/or 
cause  its  intermediate  output  to  be  performed.  Upon  completion,  the  same 
routing  message  is  typed  on  the  screen.  Entering  0 to  return  sets  array 
XKB  equal  to  the  work  array  XC  and  causes  a return  to  the  main  program. 

Entering  1 to  plot  causes  the  following  to  happen.  Subroutine 
MAG  takes  the  complex  value  of  the  input  compensator  XKB  stored  in  array 
XC  and  subroutine  MCMP  performs  the  product  of  these  complex  values  with 
those  values  in  ^ at  the  entrance  of  the  routine,  stored  in  COMMON  /A33/. 

This  matrix  product  is  stored  in  the  Q arrays  in  COMMON  /A4/.  Subroutine 
GERSH  then  computes  the  moduli  of  the  compensated  plant  and  subroutine  DRW 
is  called  at  entry  point  DRWl  for  a graphical  display  of  the  computed 
information.  After  all  the  desired  plots  are  inspected,  the  original  routing 
message  appears  on  the  screen  allowing  for  further  modification  or  a 
return  to  the  main  program. 

The  effects  of  repeative  changes  in  array  XC  are  not  accumulative 
and  those  values  of  the  system  stored  in  the  Q arrays  in  COMMON  /A4/  are 
for  the  last  input  compensator  used.  Caution  should  again  be  exercised 
when  leaving  the  package  and  entering  another  package  from  the  main  program. 

4.4.3.  Subroutine.  CLSP 

Typing  3 for  the  closed  loop  package  option  in  the  main  program 
calls  subroutine  CLSP.  This  programming  segment  allows  for  the  interactive 
modification  of  the  feedback  matrix,  F,  the  computation  of  the  Ostrowski 
bands,  and  the  graphical  display  of  the  resulting  closed  loop  composite 


inverse  Nyquist  diagrams  using  the  maximum  possible  screen  size.  The 
functional  flow  chart  of  this  package  is  included  as  Figure  4.4.3. 
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When  entering  the  routine,  work  matrix  XS  is  set  equal  to  the 
value  of  the  feedback  matrix  F from  COMMON  /Al/.  Subroutine  GERSH  is 
called  to  compute  the  Gershgorin  bands  and  subroutine  OS TROW  computes  the 
value  of  the  Ostrowski  band  radii  from  the  Gershgorin  bands.  The  message 

DEALING  W/  CLOSED  LOOP  SYSTEM 

is  typed  on  the  terminal  screen  to  remind  the  user  that  the  matrix  XS  has 
taken  on  the  function  of  the  feedback  matrix  F. 

The  routing  message 

TYPE:  0 TO  RETURN,  1 TO  PLOT,  2 TO  MODIFY 

appears  on  the  screen.  Entering  0 sets  matrix  F equal  to  work  matrix  XS 
and  a return  is  made  tc  the  main  program.  Entering  2 to  modify,  results  in 
the  calling  of  subroutine  MODXS  for  changing  values  in  the  feedback  matrix 
and/or  their  output.  The  return  from  routine  MODXS  causes  the  original 
routing  message  to  again  appear  on  the  terminal  screen. 

If  a graphical  output  is  desired  the  integer  1 is  entered.  This 
causes  subroutine  OSTROW  to  compute  the  Ostrowski  bands  and  subroutine 
PTBIGl  is  called  for  the  graphical  investigation  of  the  diagonal  closed 
loop  inverse  frequency  loci.  Returning  from  the  graphics  routine  causes 


the  original  routing  message  to  be  typed  on  the  screen  for  further 
modification  of  feedback  gains  or  a return  to  the  main  program  for  a different 
design  option  to  be  executed. 


This  package  is  also  useful  in  showing  the  effects  of  transducer 
failures.  By  changing  a specific  feedback  element  from  its  design  value  to 
zero  with  the  interactive  modifying  routine  MODXS,  the  composite  frequency 
loci  of  the  system  with  this  failure  may  be  investigated. 


4.5.  Discussion  of  Autopilot  Example 

The  autopilot  plant  transfer  function  (4.2.8)  and  (4.2.9)  and  the 
programming  variables  in  COMMON  /A2/  were  assigned  in  data  subroutine 
DATA12.  Constant  input  compensator  XKA  was  initialized  to  the  identity 
matrix  as  well.  The  open  loop  package  OLISP  was  used  to  investigate  the 
composite  frequency  loci  of  the  plant. 

The  plant  without  constant  matrix  compensation  is  observed  to  be 
highly  nondominant  (Figures  4.5. 1-4. 5. 7).  Several  different  constant 
input  compensator  matrices  were  constructed  in  order  to  achieve  diagonal 
dominance  of  the  open  loop  system.  None  of  these  attempts  were  successful 
in  achieving  a diagonally  dominant  open  loop  system  at  all  frequencies. 

In  fact,  following  Hankins  [13]  it  can  be  shown  that  the  auto- 
pilot plant  cannot  be  put  into  a diagonally  dominant  form  with  a constant 
matrix  compensator . 

One  attempts  to  obtain  row  diagonal  dominance,  as  this  has  certain 
advantages  when  matrix  compensation  is  placed  in  the  forward  path  [14]  and 
the  system  was  originally  constructed  so  this  would  be  the  case.  Therefore 

A 

K is  partitioned  by  rows  and  each  of  these  rows  is  designed  by  considering 

A 

the  row  dominance  in  that  corresponding  row  of  Q.  This  requires  G(jou)  to 


PLOT  BMNC»S?  TYPE=  O FOR  NO.  1 FOR  ROWS. 2 FOR  COLUMNS 


with  Gershgorin  bands  for  0<(»)<1.6. 


Figure  4.5.7.  §22  with  Gershgorin  bands  for  0<u)<4.0. 


77 


I 

f 

[ 

i 


( 


I- 

I 

I 

t 

I 

i 

t 


c 

i 


I 

I 

I 

I 

I 

I 


be  partitioaed  by  columns  as 

G(jcu)  = [o',^+je^:Q'2+jP2^'  (^.5.1) 

Hawkins  shows  that  diagonal  dominance  cannot  be  obtained  if  the 
two  eigenvalues  of  the  matrix 

C,  = fa.aT+P.pT  - a.aT  - P.pT]  (4.5.2) 

i 11  11  j j j j V 

have  the  same  sign.  He  then  carries  this  idea  further  to  a generalized 
graphical  method  of  obtained  dominance  for  a two  input  two  output  plant 
by  specifying  numerical  bounds  for  the  elements  in  such  a constant  compen- 
sating matrix. 

For  the  autopilot  plant,  the  value  of  the  inverse  transfer 
function  at  ju)=  .1  is 


G(j.l) 


-(2.40+J1.95)  -(1.48  + jl.73) 

(.74-j.94)  (.56-j.66) 


(4.4.3) 


Although  is  already  both  row  dominant  and  column  dominant 
for  this  frequency,  is  computed  from  (4.5.2)  as 


C 


1 


3,18 

-1,05 


-1,05 

,608 


(4.4.4) 


Observe  that  is  symmetric  by  construction  and  that  C2  = -0^^.  The  eigen- 
values of  are  found  to  be  +3.47  and  +.32.  As  they  have  the  same  sign, 
diagonal  dominance  in  this  manner  can  never  be  achieved  by  a constant 

A 

compensator  K. 
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As  open  loop  diagonal  dominance  cannot  be  obtained  with  constant 
compensation,  a design  for  the  autopilot  regulator  problem  could  proceed 
in  one  of  the  following  directions: 

1)  Nothing  prevents  one  from  attempting  to  construct  a nondiagonal 
input  matrix  compensator,  which  is  a function  of  frequency,  to  obtain  an 
open  loop  diagonally  dominant  system.  However,  one  should  keep  in  mind  the 
necessary  complexity  of  the  control  scheme  which  will  result  after  the 
additional  diagonal  compensators  for  the  desired  input-output  time  domain 
responses  are  included.  Present  controlling  methods  in  actual  use  should  be 
surveyed  to  see  if  these  means  are  justifiable  in  terms  of  maintenance 

and  realiability . The  resulting  improvement  in  performance  of  the  operating 
system  as  a function  of  control  will  ultimately  be  the 

deciding  factor.  Clearly  more  fundamental  work  in  achieving  open  loop 
dominance  via  these  means  is  needed. 

2)  Off  diagonal  elements  which  are  a function  of  frequency  may  be 
inserted  in  the  feedback  path  to  allow  for  transient  compensation  of  the 
interaction  of  the  transfer  functions  in  the  closed  loop  system.  As  the 
inverse  Nyquist  method  has  been  used  only  in  industrial  control  problems 
which  have  not  needed  this  device,  all  the  consequences  of  this  type  of 
control  are  not  fully  known. 

3)  We  chose  as  outputs  from  the  state  description  of  the  plant  the 
yaw  angle  and  the  roll  angle.  These  motions  are  highly  coupled  through 

the  dynamics  of  the  aircraft  with  the  sideslip  angle.  An  analysis  should  be 
performed  to  see  how  the  sideslip  angle  could  be  Included  as  an  output  so 
that  the  resulting  description  of  the  inverse  plant  has  a dominant  form. 


Unfortunately,  an  algorithm  for  such  an  analysis  has  not  yet  been  determined, 
4)  It  should  be  kept  in  mind  that  one  was  given  highly  accurate  but 
typical  values  in  the  perturbed  stated  equations  about  an  operating  point 
for  the  plant  description.  These  resulted  in  frequency  loci  which  were 
somewhat  atypical  in  that  they  exhibited  peculiar  behavior  in  certain  narrow 
frequency  regions.  The  validity  of  this  precise  energy  description  of  the 
plant  should  be  questioned.  It  of  course,  would  be  better  to  measure  the 
frequency  description  of  the  actual  plant  in  operation.  This  could  be  done 
at  various  altitudes  and  velocities,  and  under  different  loading  conditions. 
Perhaps  a more  meaningful  place  to  begin  a design  would  be  an  investigation 
of  the  change  in  the  input -output  frequency  description  under  such  different 
operating  conditions. 
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5.  CONCLUDING  REMARKS 

The  inverse  Nyquist  array  method  represents  an  extremely  practical 
way  in  which  one  may  solve  multivariable  control  design  problems.  It  is 
somewhat  limited  in  that  the  number  of  system  inputs  and  outputs  must  be 
equal,  A large  number  of  plants  fit  into  this  category,  however. 

Certainly  an  important  aspect  of  the  method  is  the  great  freedom  the 
designer  has  in  formulating  his  problem.  This  is  r.ot  only  from  a standpoint 
of  the  structure  of  the  controlling  scheme.  System  specifications  in  terms 
of  marginally  acceptable  performance  may  be  included  in  the  event  of  controller 
variations  on  total  failure  or  transducer  and  actuator  hardware.  Plant 
variations  may  be  dealt  with  in  a straightforward  manner. 

The  results  presented  here  deal  only  '*it:>  a linear  and  time  invariant 
plant.  Some  theory  for  the  nonlinear  and  time  varying  case  lased  on  describing 
functions  does  exist  [2,  pp,  179-184],  These  results  could  probably  be  extended 
to  sample  data  systems  as  well. 

The  somewhat  intuitive  action  of  the  design  method  naturally  calls 
for  an  interactive  computer  with  graphics  factiity  for  its  solution.  The  race _t 
development  and  present  availability  of  this  hardware  makes  such  an  alternative 
m.ethod  possible.  In  addition,  it  sheds  a new  light  cn  graphically  oriented 
design  procedures  in  general. 

This  method  has  been  applied  to  a variety  of  practical  applications 
[2,  pp.  190-219],  As  it  becomes  more  widely  accepted,  results  will  come  forth 
allowing  for  additional  theoretical  work.  In  fact,  very  recent  developments 


allow  the  diagonal  dominance  criterion  to  be  considerably  loosened;  the  need 
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for  detailed  high  frequency  performance  in  the  inverse  plant  has  been  elimi- 
nated as  well  [14].  Without  doubt,  further  work  using  this  alternative  type 
of  approach  will  yield  practical  solutions  to  many  complex  design  problems 
in  the  near  future. 
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APPENDIX  1 


The  data  file  Ql.DAT  has  a zero  in  the  second  space  of  the  first 


line  and  may  be  made  using  the  CREATE  command  as  follows: 


CREATE  Ql.DAT 

Create:  DSKC:Q1.DAT[404,443] 
00100  0$ 

*ES 

[DSKC:Q1.DAT[404,443]] 


TYPE  Ql.DAT 
0 


The  main  program  stored  in  disk  area  [404.443J  may  be  executated 
from  the  monitor  command 


EXB@MAIN.CMD 

where  the  command  file  MAIN.CMD  is 


AUTOMP.FOR 
.OLISP.FOR 
.IFDCP.FOR 
,CLSP.FOR 
, COMP. FOR 
.GRAPH.  FOR 
.AUX.FOR 
.DATA  12 .FOR 
.QINV2 .FOR 
.SYS :AG2 10/SEA 
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This  procedure  creates  relative  binary  files  of  'FILENAME ' .REL  as 


the  following; 


EXE@MAIN.CMD 


FORTRAN: 

MAIN. 

FORTRAN: 

OLISP 

FORTRAN: 

IFDCP 

FORTRAN: 

CLSP 

FORTRAN: 

MAP 

MAG 

MIMP 

MCMP 

CaiRSH 

OS TROW 

MMIJ 

FORTRAN : 

PTNYl 

PTGRl 

PTOSTl 

PTBIGl 

MAK4 

DRW 

FORTRAN: 

MODIA2 

MODXS 

MODXC 

QOUT 

QDSKWR 

QDSKRD 

OUT 

FINISH 

FORTRAN: 

DATA 12 

FORTRAN : 

QINV 


AUTOIP 


OLISP 


IFDCP 


GRAPH 


DATA 12 


Q1NV2 


j 
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If  these  relative  binary  files  are  already  created  and  data  is 
stored  in  disk  file  Ql.DAT,  the  executation  of  the  main  program  and  its 
termination  appears  as 


EXBCaMAIN.CMD 
LINK:  Loading 

[LNKXCT  AUTOMP  Execution] 

Q VALUES  WERE  OBTAINED  FROM  DISK 

TYPE:  0 TO  STOP 

TYPE:  1 FOR  OPEN  LOOP  PACKAGE 

TYPE:  2 FOR  INPUT  COMPENSATOR  PACKAGE 

TYPE:  3 FOR  CLOSED  LOOP  PACKAGE 

TYPE:  4 FOR  OUTPUT 

TYPE:  5 TO  MODIFY  INTIGER  COMMON /A2 /VALUES 
0 

YOU  ARE  DONE: BEST  PRINT  OUT  SYSTEM'.'.'. 

FOR  OUTPUT  TYPE; 

0 FOR  NONE,l  FOR  TTY, 2 FOR  LP 


0 

STORE  Q VALUES?  TYPE  0 FOR  NO,  1 FOR  YES 
0 

5:  @STOP 

END  OF  EXECUTION 

CPU  TIME:  3.80  ELAPSED  TIME:  1:30.37 
EXIT 


APPENDIX  2 


The  twenty-one  subroutines  are  stored  individually  on  disk  in 
area  [404,443]  in  file  names  corresponding  to  the  subroutine  name  with 
extension  of  SAV.  They  have  been  left  unprotected  and  therefore  may  be 
copies  into  a user's  disk  area  by  ececuting  the  monitor  command 

COPYXXXXXX . SAV [404 ,443] 

where  XXXXXX  is  the  appropriate  filename.  One  may  copy  all  twenty  one 
routines  at  once  by  using  the  option 

COPY*. SAV [404,443] 

In  a similiar  fashion  the  programming  for  the  autopilot  problem 
may  be  transferred  to  a different  users  area.  These  files  include: 

I.  COMP. FOR  which  contains  the  computational  routines: 

1)  Subroutine  MAP,  which  computes  points  on  the  frequency  loci  for 
any  region  on  the  Nyquist  contour  D. 

2)  Subroutine  MAG,  which  computes  the  complex  value  of  any  transfer 
function  matrix. 

3)  Subroutine  MIMP,  which  takes  the  product  of  a real  matrix  and  a 
complex  matrix  stored  by  real  and  imaginary  parts. 

4)  Subroutine  MCMP,  which  takes  the  product  of  two  complex  matrices 
each  stored  by  real  and  imaginary  parts. 


5)  Subroutine  GERSH,  which  computes  various  sums  of  the  moduli  of 


elements  in  the  open  loop  plant. 
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6)  Subroutine  OSTROW,  which  computes  the  radi  of  the  Ostrowski  bands 
in  the  closed  loop  system, 

7)  Subroutine  MMIJ,  which  searches  for  minimum  and  maximum  values  on 
a composite  frequency  locus •' 

See  Table  A. 2.1  for  a listing. 


I 

I 

I 

I 

I 

I 

I 

I 


II.  GRAPH, FOR  which  contains  the  graphics  routines: 

1)  Subroutine  PTNYl,  which  plots  a single  frequency  loci  on  the  screen. 

2)  Subroutine  PTGRl,  which  plots  the  Gershgorin  bands  based  on  rows 
or  columns  on  the  screen. 

3)  Subroutine  PTOSTl,  which  plots  the  Ostrowski  bands  based  on  rows 
or  columns  on  the  screen. 

4)  Subroutine  PTBIGl,  which  interactively  causes  a single  composite 
Nyquist  diagram  with  Ostrowski  bands  using  the  maximum  possible 
screen  size  to  be  plotted. 

5)  Subroutine  MAK4,  which  plots  four  separate  frequency  loci  with 
Gershgorin  bands  on  the  screen  simultaneously, 

A listing  of  this  program  is  given  in  Table  A. 2. 2. 

III.  AUX.FOR  which  contains  the  auxiliary  routines: 

1)  Subroutine  M0DIA2,  which  interactively  modifies  the  range  of  the 
frequency  loci  and  the  number  of  Gershgorin  or  Ostrowski  bands. 

2)  Subroutine  MODXS,  which  interactively  modifies  and/or  outputs 
scalar  matricies. 

3)  Subroutine  MODXC  interactively  modifies  and/or  outputs  frequency 
dependent  compensators. 


Table  A. 2.1.  COMP.  FOR  listing 
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I 

I 

I 


e^ll9a 

002ae 

00300 

0oa0o 

00S00 

00600 

00700 

00600 

00900 

01000 

01100 

01200 

01300 

01400 

01S00 

01600 

01700 

01600 

01900 

02000 

02100 

02200 

02300 

02400 

02500 

02600 

02700 

02600 

02900 

03000 

03100 

03200 

03300 

03400 

03500 

03600 

03700 

03800 

03900 

04000 

04100 

04200 

04300 

04400 

04500 

04600 

04700 

04600 

04900 

05000 

05100 

05200 

05300 

05400 

05500 

05600 

05700 

05600 

05900 

06000 


Su9»00TIN£  MAPdCIRCL) 


C 

970 


720 

730 


740 


750 

751 
C 


TAKES  COMPUX  HAP  or  G InTq  ()(«l,Xlt) 

ICIRCL  • 0 FOR  S BETWEEN  ONEGO  ANq  OMEGND  ONUY 
■ I FOR  S ON  THE  SEMICIRCULAR  ARC  ONLY 
• a OR  ANY  OTHER  value  FOR  ENTIRE  NYOUIST  PATH 

implicit  COmPlExCCI 

COMMON/A2/nMe&0,OMEGNO,OTOMEG,ISTART, IENO,NP,MO 


USES  COMMo^/il /,C0MM0N/a32/  ALSO 


742 


COrtMON/Al/GC^.^*t'*^)  ,XL6C2,2,3,2), 

1 XKA(2,2),xUA(2,2),F(2,2) 

COMMON/ A 32 /XRl (2,2,2011 , XII (2,2,201) 

• «*K**«>*******i>»**»*(***r>***»***********»**********l*l 

FDRMAT(X,'*  pole  at  omega  » ',F7,2,2X, 'FOR  G(',I1,',', 
II. ',',13,')') 

NTIMES»  ( (OMEGNO-OMEG0) /OTOMEC)^! 

00  760  I»1,M0 
00  780  J»1,M0 
IF(ICIRCL.EQ.I)  go  to  751 
Maps  top  strip  of  jw  axis 
Do  750  N=1,NTIME5 
IF(N.GT.IEND)  go  to  780 
OmEGA«OMEG0+(N-1) aOTOMEC 
SIGMAI0.0 
GO  TO  730 
SlGMAt.OTQMEC 
type  970, omega, I, J,N 
css »CmPlx  (SIGMA, OMEGA) 

CNUM* (0,0,0,0) 

COOMP(0, 0,0,0) 

DO  740  K»2,6 

CNUMiCNUM+CmPlX (G(I, J,K, 1) ,O.0)*C5Sa»(K-1) 
CUOMicnOM*CMPLX(G(I, J,K,2) ,0.0)*CSS»»(K-1) 

CONTINUE 

CNI)M.cnUM*CMPLX  (G(I,  J,  1 , 1)  ,0.0) 

.Cf)OMicDOM*CMPLX(G(I,  J,l,2),0.0) 

IF(CDOM,EO.(0,0,0,O))  GO  TO  720 

CpCNUM/COOM 

XRl  (I,J,N)»rEAL(C) 

XII  (I,J,n)«AIMAG(C) 
continue 

IF(ICIRCL.Eq,0)  go  to  780 
IFCICIRCL.ElJ.l)  NTIMES«0 
MAPS  SEMICIRCULAR  RADIUS  PART 
00  760  N»l,99 
N1»N*NTIMES 

OmeGA.PmEGND»COS(Na,0314159) 

SIGmAiOmEGNo»SIN(N«,0314159) 

CSSpCmPLX (SIGMA, OMEGA) 

Cnump (0,0, 0,0} 

COOMl (0.0,0,0) 

00  V42  K»2,6 

CNUMaCNuM»CMPLX(G(l, J,K, 1) ,0,0)»CSS»*(K-D 
CDOM.cOOMaCMPlX (G(I, J,K,2) ,0.0) aCSSaaCK-I) 

continue 


I 
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I 


B6100 

06200 

06300 

06000 

06500 

06600 

06700 

06900 

760 

06900 

07000 

C 

07100 

07200 

07300 

07900 

07500 

07600 

723 

07700 

07900 

733 

07900 

09000 

08100 

09200 

09300 

09400 

743 

09500 

09600 

09700 

09900 

09900 

09000 

09100 

770 

09200 

790 

09300 

09400 

09500 

C 

09600 

c 

09700 

c 

09900 

c 

09900 

c 

10000 

i' 

10100 

10200 

10300 

c 

10400 

c 

10500 

c 

10600 

10700 

10900 

c 

10900 

c 

11000 

c 

11100 

11200 

11300 

c 

11400 

11500 

11600 

11700 

11900 

11900 

12000 

20 

CNUM.CNUM  + CMPLX  CG(I, J, 1 , 1) ,0,0) 

C00M»C00«  + CMPLX  CR  Cl , J, 1 , 2) ,0.0) 
IFtCOOM.Ea.(0,0,0,0))TYPE  970,OMEG*,I» J»N 
IFCCOOM.EO,  (H, 0,0,0))  GO  TO  760 
CtCNUM/COOM 

XRl  Cl,  J,^*fJTIMES)»RE*LCC) 

Xll  Cl,  J,*"»NTIMES)«*IM*GCC) 

continue 

IFCICIRCL.EQ.1)  CO  TO  790 

M*ps  bottom  strip  of  jw  axis 

Do  770  n*1,nTIM£S 

N2»N+99+NTIHES 

OmEGAsOmEGND-CN-1)»OTOHEG 

SIGMA*0,0 

GO  TO  735 

SICMAi.DTOMEG 

TYPE  R70,OMEGA, I, J,N2 

CSS*CMPLX (SIGHA,OHEGA) 

CNUM«CO, 0,0,0) 

CDOMi(0. 0,0,0) 

DO  7A3  K»2,6 

CNUMbCNUM  + CMPLX  CGCI,J,K,t),0,0)*CSS**CK-l) 
COOMtCOON  + CMPLX  CGCI , J,K,2) ,0,0)*CSS*«CX"i) 
CONTINUE 

CNUM«CNUM»CMPUX  CGCI, J, I , 1) ,0.0) 
CDOM«CDOM*CMPlX  CG  Cl , J, 1 ,2) ,0.0) 
IFCCDOH.EQ.C0, 0,0.0))  GO  TO  723 

C»CNUM/COOM 

XRl Cl, J,N+Rq»NTIMES)»REALCC) 

XII Ct,J,NY9qYNTlMES)»AIHAG(C) 

CONTINUE 

CONTINUE 

return 

END 


no  interactive  usage 


NO  additional  routines  called 


subroutine  mag 

TAXES  COMPLEX  value  of  XC  INTO  CXRl.XIl) 


implicit  COmplEXCC) 

common/ A? /OMEG0,OMeGNO,OTOMEG, ISTaRT, IENO, ID, MO 


USES  COmmon/a32/,COMMON/A7/  ALSO 


common/ A32/XP1 C2,2,20l) ,xn  C2,2,20l) 
COMmon/a7/XS(2,2) ,XCC2,2,6,2) 


DO  60  I*l,MO 

DO  60  J«l,MO 

DO  50  N.i.IEND 

OMEGAaOMEG0*CN-1) aDTOMEC 

SIGMAB0.0 

GO  TO  30 

SICMAi-OTOMEG 


h. 


I 

I 


12100 

12200 

900 

12300 

12900 

30 

12500 

12600 

12700 

12800 

12900 

13000 

40 

13100 

13200 

13300 

13900 

13500 

13600 

13700 

50 

13800 

60 

13900 

14000 

14100 

C 

14200 

C 

14300 

C 

19400 

C 

14500 

c 

14600 

c 

14700 

c 

14800 

14900 

c 

15000 

c 

15100 

c 

15200 

c 

15300 

15900 

15500 

c 

15600 

c 

15700 

c 

15800 

15900 

16000 

16100 


16.200 

C 

16300 

16400 

16500 

16600 

16700 

10 

16800 

16900 

17000 

17100 

17200 

17300 

23 

17400 

17500 

17600 

30 

17700 

40 

17800 

17900 

18000 

c 

Type  900,oMec*, I, j,N 

F0«maT(x,'*  pole  at  OmEC*  • »,F7,2,2X,'fOH  FUNCTION  (»,I1,',», 

11.', M3,')*) 

CSSiCHPl*  (SIGMA, OMEGA) 

Cnum»(0. 0,0,0) 

cnnM« (0,0, 0,0) 

DO  ao  Ki0,b 

CNUM»CNUm*CmPl*(*C(I, J,K, 1) ,0,0)»CSS*»(K-l) 

CnOMacUOMtCMPuA (XC (I, J,K,2) ,0,0)*CSS«»(K"1) 

CONTINUE 

CNUH.CNuM*CmPLX(XC(I, J,l, 1) ,0,0) 

CUOM*COOM*CMPLX(XC(I, J, 1,2),0,0) 

IF (COOH.CQ, (0,0,0,0))  GO  TO  20 

C»CNUM/COOM 

XPl (I, J,N)«SEAL(C) 

xII(I,J,N)«AIMAG(C) 

CONTINUE 

Continue 

RETURN 

END 


NO  interactive  usage 


NO  additional  routines  called 


Subroutine  mimp 


TAKES  matrix  MULTP,  of  XSaCMPLX (XR1,XI1) 
STORES  the  result  in  0 


implicit  COmPlEX(C) 

common/ A 2/0 1 ,02,03, 1 ST ART, I END, 10, MO 


USES  COmmon/A5I/,COMMON/A32/,COMMON/A«/,COMMON/AT/  ALSO 


C0MM0N/A3WC1  (2,2)  ,C2(2,2)  ,S1  (2)  ,S2(2) 
common/ A32/XR I (2,2,201) , XI 1 (2,2,201) 
C0MM0N/A4/0R (2, 2, 201), 01 (2,2,201) 
COMmON/A7/XS(2,2) ,XC(2,2,6,2) 


00  90  N»l,I£NO 
DO  10  I»1,M0 
00  10  J»1,M0 

Cl (I, J)«CMPlx(XR1 (I,J,N) ,XI1(I,J,N)) 

continue 
Oo  30  I«1,M0 
do  30  J«1,M0 
Cx«(0, 0,0,0) 

00  20  K«1,M0 
CX«CX«XS(1,K)*C1(K,J) 

CONTINUE 

QRCI,J,N).REAl(CX) 

QI (I, J,N).AXMaG(CX) 

CONTINUE 

Continue 

RETURN 

ENO 


18100 

C 

18200 

C 

18300 

C 

1BR00 

C 

18500 

C 

18600 

18TO0 

C 

18800 

c 

18000 

c 

19000 

c 

19100 

19200 

19300 

c 

19R00 

c 

195C0 

c 

19600 

19700 

19800 

19900 

20000 

c 

20100 

20200 

20300 

20400 

20500 

20600 

c 

20700 

20800 

20900 

21000 

20 

21100 

21200 

21300 

30 

21400 

40 

21500 

21600 

21700 

50 

21800 

C 

21900 

C 

22000 

C 

22100 

c 

22200 

c 

22300 

c 

22400 

22500 

c 

22600 

c 

22700 

c 

22800 

c 

22900 

23000 

23100 

c 

23200 

c 

23300 

c 

23400 

23500 

23600 

23700 

c 

23800 

23900 

24000 

c 

NO  interactive  usage 


NO  AbOITlONAL  ROUTINES  CALLED 


subroutine  mcmp 


TAKES  matrix  MUlTP,  of  CMPLX(X»l,Xin»CHPLX(XR2,XI2) 
STORES  result  In  a 


IMPLICIT  COmPlEx(C) 

COMMON/A2/ni,oa,n3,ISTART,IENO,ID,MO 


USES  C0MMQN/A3J /,COMMON/AJ2/,COMMON/A33/,COMMON/Aa/  ALSO 


common/ a3i /Cl (2,?) ,cac2,2) ,sj  tai isatz) 

COMMON/ A3a/XR1 (2,a,20n , XII  (2,2»20l) 

coMMON/A33/xwata,a,a0n ,xiaca, 2,201) 

COMMON/Aa/OR(2,2,201) ,01 (2,2,201) 


DO  50  N*I,IEND 
DO  20  1*1, MO 
DO  20  J«1,M0 

Cl (I, J)*CMPlX(XR1 (I, J,N) ,xn (I, J,N)) 
C2(I,  J)»CMPi.X(XR2(I,  J,N)  ,XT2(I,  J,N)) 
CONTINUE 
Cx*(0, 0,0,0) 

DO  30  K»1,M0 
CX*CX4CI(I,K)*C2(K,J) 

CONTINUE 

0R(I,J,N)*»EAL(CX) 

Oin,J,N)*AlMAC(CX) 

CONTINUE 

continue 

return 

END 


NO  interactive  usage 


NO  additional  routines  used 


subroutine  gersh 


computes  GERSHCOHIN  bands  for  current  q values 
PLACES  result  in  C0MM0N/A5/  BANOS 


implicit  CDmPlEXIC) 

COMMON/A2/OI,D2,03,I8TART,IENO,ID,MO 


USES  C0MM0N/*3l/,C0MM0N/Aa/,L0MM0N/A5/  ALSO 


COMMON/ A3 1 /Cl (2,2) ,C2(2,2) ,S1 (2) ,S2(2) 
COMMON/A«/QR(2,2,20l) ,01  (2,2,201) 
COMMON/ A5/B AND (2, 2, 20 1 ,2) 


DO  30  N*1,IEN0 
DO  10  1*1, MO 
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1 

I 


241B0 

24200 

24300  10 

24400  * 

24500 

24600 

24700 

24000 

24900 

25000 

25100  20 

25200 

25300 

25400  30 

25500 
25600 
25700  C 
25000  C 
25900  C 
26000  C 
26100  C 

26200  r 

26300  C 
26400 
26500  C 
26600 
26700  C 
26000  C 
26900 
2T000 
27100  C 
27200  C 
27300  C 
27400 
27500 
27600 
27700 
27000  C 
27900 
28000 
20100 
20200 
20300 
20400 
20500 
20600 
28700 
20600 
20900 
29000 

29100  20 

29200 

29S00 

29400  30 

29500  40 

29600 
29700 
29000  C 
29900  C 
30000  C 


00  10  J«1|M0 

can,  j)iCMPui(0R(i,  jfN) 

continue 

00  30  I*1,M0 
00  30  J«1,H0 
XR«O,0 
XC«0.0 

DO  20  K«1,M0 

IFtK.NE.J)  XR»XR+C*BS(C2(I,K)) 
IF(K.NE.I)  XC«XC*C*B5(C2(K, J)) 
CONTINUE 

BANDCI, J,N,1)«XR 
BANUCI, J,N,2)*XC 

continue 

return 

f,NO 


NO  interactive  usage 


NO  additional  routines  used 


subroutine  ostrow 


COMP.  OSTROmSKI  bands  from  current  GERSH,  BANDS, 0 AND  F VALUES 
places  result  in  C0MM0N/A6/  OST 


implicit  COmPlExCC) 

COMMON/a2/OMEG0,OMEGND,DTOMEG,ISTART, IENO,NP,MO 


USES  COMMoN/A7/,COMMON/A«/,COMMON/A5/,COMMON/A6/  also 


COMMON/A«/OR(e,2,20n,OI(2,2,201) 
COMMON/A5/BAN0  12,2,201, 2) 

COMMON/A6/OST(2,P01,2) 

C0MM0N/A7/XS(2,2),XC(2,2,6,2) 


00  40  Nil,IENO 

XORa0,0 

XOC*0.0 

00  30  III, MO 

00  20  Jal,HO 

IFCJ.EO.U  GO  TO  20 

C>CMPLX((XS(J, J}40R(J, J,N}),OX(J, J,N)) 
HaCABSIC) 

X0R1«BAnD(J, J,N,1)/H 

XOCl«BAND(J, J,N,2)/M 

IF(XOHl.GT.XOR)  X0R»X0R1 

IF(XOCI.GT.XOC)  XOCiXOCl 

CONTINUE 

OST(I,N, IjaxOR 

0ST(I,N,2)aX0C 

CONTINUE 

CONTINUE 

RETURN 

END 


NO  interactive  usage 


V 


NO  additional  HQUTINES  USED 


30100  C 

30000  C 

30300  C 

30400  C 

30S00 

30fr00  C 

30700  C 

3OS00  C 

30900 

3(000  C 

31100  C 

31200  C 

31300 

31400 

31500  C 

3lb00 

31700 

31S00 

31900 

32000 

32100 

32200 

32300 

32900 

32300 

32600 

32700 

32S00 

32900 

33000 

33(00  20 

33200 
33300 
33400  C 
33500  C 
33600  C 
33700  C 
33600  C 
33900  C 


Subroutine  mmijci,j,xhin,xmax,vmin,tmax) 


SEAhCHS  Cl.J)  COMPOSIT  NYOUIST  DIAGRAM  FOR  (X,T)  MIN/MAX 


COMMON/ A2/01 ,02,03, 1ST 4RT, I END, NIDI, MO 


USES  common/a4,common/a5/  also 


COMMON/ A4/0R (2, 2, 2011 ,01(2,2,201) 
common /AS/ BAND  (2, 2, 201 ,2) 


DO  20  NL41,IEND 
1RC«1 

1F(BAND(I, J,NL, 1) ,GE,BAND(I, J,NL>2))  IRCa2 
XMINl.QRd,  J,NL) -BAND!  1,J,NL,  IRC) 
XMAxl«OR(I, J,NL)*HANO(I, J,NL,1RC) 
IF(I.NE.J)  XMINi.qr(I, j,nl) 

If(i.ne.J)  xMAXi«afi(i,j,NL) 
IFtxMINl.UT.XMlN)  XMINiXMINl 
IF(XM*XI,GT.XMAX)  xmax*xmaxi 
TMINl.QI (I, J,NL)-BAN0(I, J,NL,IRC) 
YMAX1«QICI, J,NL)4BAND(I, J,NL,IRC) 
IFCI.NE.J)  YMINl.Qt  (I, J,NL) 

IF(I.NE.J)  YM*X1«QICI,J,NL) 

IFCyMINi.lt, ymIn)  ymIN«yMIN( 

IF (YM*X1 .GT.YMAX)  YMAX»Y«AX1 

continue 

RETURN 

ENO 


NO  interactive  USAGE 


NO  additional  routines  called 
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1 

t 
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Table  A. 2. 2.  GRAPH. FOR  listing 


(■ 

t 

I 

i 

f 

t 

i 

i 

I 

( 


22120 

C 

22200 

C 

02322 

22422 

C 

20522 

C 

00622 

C 

00T20 

00800 

C 

20422 

c 

01222 

c 

21100 

01200 

c 

01320 

01420 

c 

01500 

c 

01600 

01T20 

21820 

21402 

02020 

02102 

02200 

10 

02300 

02400 

02500 

02600 

02T00 

c 

02800 

c 

22420 

c 

03000 

c 

03100 

c 

03200 

c 

03300 

c 

23400 

03520 

c 

03600 

c 

03720 

c 

03820 

c 

03400 

c 

04200 

24102 

c 

04200 

c 

04300 

04400 

24500 

c 

04600 

04700 

04802 

04400 

25200 

05100 

05200 

05300 

05402 

05500 

to 

25620 

25700 

20 

05802 

25400 

06020 

c 

Su9HOuTI^t  PTNYJ (1. J) 


Plots  (I,jt  nyquist  oi*c«*m 


co«''ONi/Aa/ni,02,03,isr*PT,ieNo,io,MO 


USES  CO»»QH/^a/  ALSO 


co«Mou/*tt/ORca.2,20n  ,01  ta.s.apn 


dime^^sio*-  Ync202),xica0a5 

note  I FOW  PU0TTIS5  f*0«E  THAN  221  POINTS  IN  Q CHANGE 
OIhEnSIONS  ACCOHDInCLY 
lMANY«TENO-ISTAaT»l 
XS (1) ilMANY 
XlCn*IMANY 
00  12  N»1,IMANY 
XR(N.n»0PtI,J,^*I9TABT-l) 

XI  (N.n  ■Old,  J,N.IST*HT-l) 

CONTINUE 

CALL  CHECK(xR,xn 
CALL  DSPLAY(XR,XI) 
return 
end 


NO  interactive  USAGE 


NO  additional  routines  used 


Subroutine  ptgri (i , j, irc) 


PLOTS  GERSHGORiN  PANDS  OF  (I.J)  ENTRY  OF  0 

IRC  » 1 FOR  Plotting  bands  based  on  robs 
■ 2 FOR  plotting  bands  based  ON  columns 


C0HH0N/A2/0mEC2,DhEGND,DT0HEG, ISTART, IEND, NP, MO 


USES  COMM0N/Aa/,COMMQN/A5/ 
C0HM0N/A«/tlR(2,2,2?n  ,QI  (2,2,221) 
COMMON/ AS/ 8 AND (2, 2, 22 1 ,2) 


DIMENSION  XC102) ,Y(122) 

NSTOP*(IENO-ISTART*n/NP*l 

X(t)*I21 

YC1)*I01 

00  20  n.i,nsTOP 

NI» (N-i ) .NP.ISTART 

DO  10  F>1,101 

X(K«l)«OR(I,J,Nn*BAND(I,J,NI,IRC)*COS((K>n*,0b283l6) 
Y(aa1}iQI (I, J,NI)*aANO(I, J,NI,IRC) ASlN((K.n*,2b2SS18) 
CONTINUE 
CALL  CPLOT(X,Y) 

CONTINUE 

RETURN 

ENO 


L. 


J 


d 


1 

E 
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Obiee  c 

06200  C 
06300  C 
06000  C 
0650U  C 
06600  C 
06T00 
06800  C 
06900  C 
07000  C 
0Tt00  C 
07200  C 
07300 
07000  C 
07500  C 
07600  C 
07700 
07800 
07900  C 
08000 
06100 
08200 
08300 
08000 
06500 
08600 
08700 
08800 

08900  10 

09000 

09100  20 

09200 
09300 
09000  C 
09500  C 
09600  C 
09700  C 
09800  C 
09900  C 
10000  C 
10100 
10200  C 
10300  C 
10000  C 
10500  C 
10600  C 
10700  C 
10800  901 

10900  10 

11000 
11100 
11200 

11300  900 

11000 
11500 
11600 
11700 
11800 
11900 
12000 


NO  interactive  usage 


NO  additional  routines  called 


Subroutine  ptosti ci, ircj 


Plots  dstrowski  bands  for  h(i,d 

IRC  C I FOR  BANDS  BASED  ON  ROWS 

■ 2 FDR  bands  based  ON  COLUMNS 


COMMnN/A2/OMEr,0,OMEGND,OTQMEG,  I ST  aR  T , IENO  , NP , MO 


USES  COMM0N/Aa/,C0MMQN/A6/  ALSO 


COMMON/A4/OR(2,2,201] ,QI (2,2,201) 
COMMON/a6/0ST(2,2O1,2) 


dimension  *(t02),YCl02) 

NST0P.(1EnD-ISTaRT*1)/NPA1 

X(l)»l0l 

r(t)«i0i 

DO  20  N»l,NST0P 

NI«(n-1)»NP+ISTART 

DO  10  Kit, 101 

X(K+l)iQR(I, I,NI)»OST(I,NI,IRC)»COS((K-1)i,0628318) 

V(K+1)«OI(I,I,NI)a.oST(1,NI,IRC)»S1N((K-1)»,0628318) 

CONTINUE 

CALL  CPL0T(X,T) 

CONTINUE 

RETURN 

END  V 


NO  interactive  usage 


NO  additional  routines  used 


Subroutine  ptbi&i 


PLOTS  NYQuIST  diagrams  W/OSTROWSKl  BANOS 
USES  maximum  SCREEN  SIZE 


NO  COMMON  BLOCK  usage 


F0RMAT(I1,X,I1) 

CONTINUE 
CALL  ANMOOE 
CALL  HOME 
TYPE  900 

FORMAT(X, 'TYPE  0,OR  TYPE  1,1/2  POR  H(l,l)  AND  OST  R0WS/C0LUM8*) 

ACCEPT  901, I, IRC 

IF(I.EO.0)  GO  TO  20 

CALL  BINITT 

CALL  NEWPAG 

CALL  SLIMX(75,1000) 

CALL  SL1MY(50,790J 
CALL  PTNYKI,!) 
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I 

r 

[ 

i 

I 

f 

f 

( 

? 


t 

I 

I 

I 

I 

I 

I 

I 


12100 

12200 

12300  20 

12«00 
12S00 
12000 
12f00 
12000  C 
12000  C 
13000  C 
13100  C 
13200  C 
13300  C 
13400  C 
13500  C 
13000  C 
13700 
13600  C 
13400  C 
14000  C 
14130 
14200  c 
14300  C 
14400  C 
14500 
14000 
14700 
14600 
14900 
15000 
15100 
15200 
15300 

15400  10 

15500 
15000 
15700 
15600 
15900 
10000 
16100 
16200 
10300 
10400 
10500 
10600 
10700 
10600 
10900 
17000  C 
17100  C 
17200  C 
17300  C 
17400  C 
17500  C 
17000  C 
17700 
17600  C 
17900  C 
16000  C 


CALL  PTOSTl  (IfIRC) 

GO  TO  10 

CONTINUE 

CALL  aNMOOE 

CALC  NEWPAG 

BETUBN 

END 


INTEBACTlVELT  ACCEPTS! 

1)  diagonal  element 

2)  0ANO  BASED  ON  RON/COLUMN 


additional  routines  called  I PTNYl.PTOSTl 


subroutine  MAK4(IKC) 


PLOTS  ALL  NyOulST  DIAGRAMS  FOR  (2X2)  SYSTEM  ON  SCREEN 


common/ A2 /Oi ,02,03 • 1D0, ID  1 , 102, MO 


NO  additional  common  USAGE 


XMlN«10000, 

YMIN»1C000, 

XMAX«-10000, 

YMAx«-10000, 

DO  10  I*1,M0 
00  10  JaifMO 
lOal 
JO*J 

call  MMU(I0,J0,XM1N,XMAX,YM1N,YMAX) 
Continue 

CALL  0LU5X(XMIN,XHAX) 

CALL  DLIMY(YMIN,YMAX) 

CALL  ERASE 
CALL  PLACE(4) 

CALL  PTNYKi.l) 

IF(IPC.nE.O)  call  PTGRl(l,l,IRC) 
call  PL*CE(5) 

CALL  PTNY1(1,2) 

CALL  PLACE(0) 

Call  ptnyk2,i) 

CALL  PLACEtT) 

CALL  PTNY1(2,2) 

1F(IRC.NE,0)  call  PTGR1(2,2,IRC) 

return 

ENO 


NO  interactive  usage 


AOOITIONAU  routines  called  t MM1J,PTNY1,PTGR1 


!l 


subroutine  ORm 
PLOT  ROuTINS  routine 
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isiee  c 
leaee  c 
ISS00  le 
1S«00 
18500 
18600 
18T00 
18800 

18900  900 

19000 

19100 

19200 

19300 

19400 

19500  901 

19600 
19700 
19800 
19900 
20000 
20100 
20200 
20300 
20400 

20500 

20600 

20700 

20800  20 

20900 

21000 

21100 

21200  902 

21300 
21400 
21500 
21600 
21700 
21800 

21900  100 

22000 
22100 
22200  C 
22300  C 
22400  C 
22500  C 
22600  C 
22700  C 
22800  C 
22900  C 
23000  C 


NO  COMMON  block  usage 


call  ANMOOE 

CALL  HOME 
CALL  NEWLIN 
CALL  NEWLIN 
Call  NEWLIN 
TYPE  900 

E0PMAT(5K, 'PLOT  OIAGBAMST  TYPEI  0 709  NO, 1 709  YES') 

ACCEPT  WflOUM 
IPnOuH.CO.O)  CO  TO  100 
ENTRY  ORWl 
CALL  ANMOOE 

Type  901 

7ormat(x, 'Enter  -0,0-  7oR  4,09  -l,J-  70r 
ACCEPT  *,I,J 

17(1. EQ.0. and, J.EO.0)  GO  TO  20 
XMINS10000O, 

YMINI100000, 

XMAXi-100000, 

YMAX*>t0000O. 

CALL  MMUn,  J,XMIN,XMAX,YMIN,YMAX) 

CALL  DLIHX(XMIN,XMAX) 

CALL  0LIMY(YMIN,YMAX) 

CALL  erase 
call  PL*CE(1) 

CALL  PTNYKIiJ) 

CONTINUE 
CALL  HOME 
CALL  NEWLIN 

type  902 

7oRmAT(5X, 'PLOT  BANOS?  TYPEl  0 709  N0,1  709  ROWS, 2 FOR  COLUMNS*) 
ACCEPT  *,IRC 

17(irc.eq.0.ano,i.ne,o)  go  to  10 
17(1. EQ.0. and, J.EO.0)  CALL  MAK4(IRC) 

17(1. E0.».*N0.J, £0,0)  GO  TO  10 
CALL  PTGRl (1, J,IRC) 

GO  TO  10 
CoNTlAliJE 
RETURN 
END 


INTERACTIVELV  ACCEPTS  I 
IJ  PLOT  YES/NO 

2)  PLOT  4 OIAGRMS/I.J  element 

3)  PLOT  GERSH. bands  nq/ROWS/COLUMNS 


AOITIONAL  routines  used  I MMIJ.PTNYl, PTGRl, MAX4 


i 


/ 
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4)  Subroutine  QOUT  outputs  the  numerical  value  of  the  open  loop  system 

A 

matrix  Q(s)  to  the  terminal  screen  or  linepr inter. 

5)  Subroutine  QDSKWR  interactively  stores  the  open  loop  system  Q(s) 
numerics  on  computer  disk  for  the  next  design  run. 

A 

6)  Subroutine  QDSKRD  automatically  reads  the  open  loop  system  Q(s) 
numerics  from  disk  so  they  do  not  have  to  be  recalculated  each  time 
the  program  is  run. 

7)  Subroutine  OUT  interactively  outputs  system  matrices  and  programming 
variable  to  the  terminal  screen  or  lineprinter. 

8)  Subroutine  FINISH  notifies  the  user  he  is  stopping,  causes  system 
outputs  and  disk  storage  to  possibly  be  performed,  and  terminates 
program  execution. 

A listing  of  this  program  is  given  in  Table  A. 2. 3. 

IV.  AUTOMP.FOR  contains  the  main  program  used  in  investigating  the  auto- 
pilot program. 

V.  OLISP.FOR  contains  the  open  loop  constant  input  compensator  design 
package . 

VI.  IFDCP.FOR  contains  the  open  loop  frequency  dependent  compensator  design 
package . 

VII.  CLSP.FOR  contains  the  closed  loop  feedback  design  package. 

VIII.  DATA12.F0R  contains  the  subroutine  defining  the  system  and  programming 
data  for  the  autopilot  problem. 

IX.  QINV2.F0R  contains  the  subroutine  which  takes  the  point  by  point  inverse 
of  the  autopilot  plant  data. 

See  Table  A. 2. 4 for  a program  listing. 


I 


Table  A. 2. 3.  AUX.FOR  Listing 


00100 

C 

00200 

00300 

c 

00400 

c 

00500 

c 

00600 

00700 

c 

00600 

c 

00900 

c 

01000 

01100 

c 

01200 

01300 

01400 

01500 

01600 

01700 

01600 

01900 

900 

02000 

c 

02100 

c 

02200 

c 

02300 

c 

02400 

c 

02500 

c 

02600 

c 

02700 

02600 

c 

02900 

c 

03000 

c 

03100 

03200 

c 

03300 

c 

03400 

c 

03500 

03600 

c 

03700 

03600 

c 

03900 

10 

04000 

04100 

04200 

04300 

901 

04400 

04500 

04600 

04700 

904 

04600 

20 

04900 

05000 

05100 

05200 

05300 

05400 

05500 

902 

05600 

903 

05700 

30 

05600 

05900 

06000 

40 

subsouti^e  modus 

MOOIFUS  INTECEW  values  in  C0MM0N/A2/ 
COMMDN/AS/OMEG0,nMEGND,DTOM£C, ISTaHT, IENO,NP,MO 
NO  OTHER  COMMON  BLOCK  USAGE 


type  R?0 

format (A, 'ENTER  ISTART,IENO,NP,OR  AND  RETURN*) 

ACCEPT  ..n.ia.u 

IFCIl.NE.O)  ISTART-It 

IFCI2.NE.0)  1EN0*I2 

IFdJ.NE.iM  NP.I3 

type  a,ISTART,IEND,NP 

RETURN 

END 


interactive  usage 

n ENTER  VALUES  OR  COMMA 


NO  OTHER  routines  USES 


Subroutine  modxs 

MODIFIES  ANO  OUTPUTS  MATRIX  X3 

•*********»**********«******A****< 

COMMON/A2/ni ,02,03, ID3, 101 , 102, MO 
USES  C0MMON/A7/  ALSO 

••A*****«*********«**************« 

C0MM0N/a7/XSC2,2) ,XCC2.2»b,2) 


CALL  ANMOOE 
TYPE  901 

FORMATIX. 'MODIFY  XS7  TTPEl  0 FOR  N0,1  FOR  YES*) 

ACCEPT  »,IDUM 
IF  (IOuM,EQ,0)  GO  TO  20 
TYPE  90a 

FORmATCX, 'TYPE!  1,J.X  IN  XS(I,J)»X*) 

ACCEPT 
XS(I,J)«X 
Co  TO  10 
type  9P2 

F0RMAT(X, 'OUTPUT  XS?  TYPEl  0 FOR  NC,1  FOR  TTY, 2 FOR  LP') 
ACCEPT  »,IOUM 
IF  (lOUM.EO.O)  CO  TO  40 
00  30  Ul.MO 
DO  30  Jil,MO 

IFdOuM.EO.n  type  903,  I,  J,XS(1,  J) 

IF(I0UM.£g,2)  PRINT  903, I, J,XS(I, J) 
FORMAT(X,'XS(',Il,',',n,*)  • *,Fr.2) 

CONTINUE 

CONTINUE 

RETURN 

ENO 


100 


Bbiee 

06300 

06300 

06400 

06S00 

06600 

06700 

06600 

06900 

07000 

07100 

07200 

07300 

07400 

07500 

07600 

07700 

07900 

07900 

08000 

08100 

08200 

08300 

08400 

08500 

08600 

08700 

08800 

08900 

09000 

09100 

09200 

09300 

09400 

09500 

09600 

09700 

09800 

09900 

10000 

10100 

10200 

10300 

10400 

10500 

10600 

10700 

10800 

10900 

11000 

11100 

11200 

11300 

11400 

11500 

11600 

11700 

11800 

11900 

12000 


INTERACTIVE  USAGE: 

1)  mooiev  xs  yes/no 

2)  ENTER  values 

3)  OUTPUT  XS  NQ/TTY/LP 


NO  aooitional  routines  used 


subroutine  mooxc 


MODIFIES  and  outputs  MATRIX  XC 


COMMON/ A2/D1 ,02,03, 100, lOl , I D2, MO 


USES  C0MM0N/A7/  ALSO 


C0MM0N/A7/XSC2,2) ,XCC2,2.6,2) 


CALL  ANMQDE 
TYPE  901 

EORMATIX, 'MODIFY  XC?  TYPE:  0 FOR  NO,  1 FOR  YES') 

ACCEPT  *,IDUM 
IF CIOuM.EQ.O)  go  to  20 
TYPE  902 

F0RMAT(X,'I,J,N,X,1/2  in  XCCI» J)»Xo(53«*N,  NUM/DOmO 

ACCEPT  »,1,J,N,X,N0 

XC(I,J,N»1,ND)»X 

GO  TO  10 

Continue 

type  903 

FORMATfX, 'OUTPUT  XC7  TYPE:  0 FOR  NO,  1 FOR  TTY,  2 FOR  tP') 
accept  »,1DUM 

IF (IOUM,FO.03  GO  TO  60 
IF  (I0UM.E3.2)  GO  TO  40 
Oa  30  I»1,M0 
00  30  J«1,M0 

type  904,xccr, j.3» d.xcci, j,2, n ,xc(i, j, 1, 1) 

F0RMATC11X,F7,?,3X,F7,P,3X,F7,2) 

type  905, I, J 

F0RMAT(X,'XC(',I1,',',11,'3  » ',29(»-*n 
TYPE  904,XC(I, J,3,2),XC(I, J,2,2),XC{I, J,  1,2) 
type  906 

formatcx,'  n 

CONTINUE 

IFC1DUM.NE.2)G0  to  60 
DO  50  I«1,M0 
DO  50  J«1,M0 

PRINT  904,XC(I, J,3, 1) ,XC(I, J,2,n ,XC(1, J, 1,1) 

PRINT  905, 1,J 

print  904,xc(I, J,3,2) ,XC(I,J,2,2) ,XC(I, J,1,2) 

PRINT  906 
CONTINUE 

continue 

RETURN 

END 

interactive  usage 

1)  modify  XC  YES/NO 


101 


I 

I 

I 

I 

I 

I 

I 


isioe  c 
c 

>2300  C 
124100  C 
>2S00  C 
12600  C 
12700  C 
12600 
12900  C 
13000  C 
13100  C 
13200  C 
13300  C 

n«00 

13300 
13600  C 
13700  C 
13600  C 
13900 
U000  C 

19100 

14200 

14300 

14400 

14S00 

>4600  10 

14700  20 

14600 
14900 
15000  C 
15100  C 
15200  C 
15300  C 
15400  C 
15500  C 
15600  C 
15700 
15600  C 
15900  C 
16000  C 
16100 
>6200  C 
16300  C 
16400  C 
16500 
16600  C 
16700 

16600  666 

16900 

17000 

17100 

17200 

17300 

17400  900 

17500 

17600 

17700 

>7600 

17900  901 

16000  30 


21  ENTER  VALUES 
31  OUTPUT  7C  NO/TTT/LP 


NO  ADDlTlON/kU  ROUTINES  USED 


SUHROuTIWE  OOUT(IOUT) 


Outputs  current  q values  in  common/a4/ 
lOUT  • 1 FOR  tty 
» 2 FOR  LP 


IMPLICIT  C0MPI.EXCC1 

common/ A2/nMEG0,()MEGND,OTOMEG, I START, IEnO,NP,MO 


USES  common/au/  also 


COMMON/A4/OR(2,2,2011 ,01 (2,2,201) 


00  20  I»1,M0 
00  20  Jal,MO 
DO  10  N*ISTaRT,IEND 

IF (lOUT.EQ.l)  type  •,I,J,N,OR(I,J,n),OICI,J,N) 

IF (IUUT.EU.2)  PRINT  • , I , J , N , OR ( I , J , N1 , 01 ( I , J, N) 

CONTINUE 

Continue 

RETURN 

ENO 


NO  interactive  usage 


NO  additional  routines  used 


SuRROUTINE  odsk,<r 

*****»*a***********«******««**tt«* 

writes  0 DATA  TO  OISK  FILE  01, DAT 


CQMM0N/A2/01 ,02,D3, ISTART, lENO.lDl ,M0 


USES  C0MMDN/A4/  ALSO 


common /A4 /OR (2,2,2011,01(2,2,201) 


TYPE  888 

F0RMAT(X, 'STORE  Q VALUES?  TYPE  0 FOR  NO,  I FOR  YES*) 
ACCEPT  *,IDO 
IF(IOO.EO,0)  GO  TO  40 
lOUM.l 

OPEN(UNIT.t .ACCESS* '5E00UT', file* *01. OAT') 

WRITE(l,900)  lOUM 

F0RmaT(X,I1) 

DO  30  n«i,IENO 
00  30  1*1, MO 
00  30  J*1,M0 

WRITE (1,901)  I, J,N,0R(I, J,N) ,QI (1, J,N) 
F0RMAT(X,Il,X,Il,X,13,X,Eia.T,X,El4,7) 

CONTINUE 


i 

5 


I 


I 


1 

1 

« 


102 


t 


I 

I 

I 

I 

I 

I 

f 

I 


18100 

18200 

40 

18300 

18400 

1BS0O 

C 

18800 

c 

18700 

c 

18800 

c 

18900 

c 

19000 

c 

19100 

c 

19200 

c 

19300 

19400 

c 

19500 

c 

19600 

c 

19700 

c 

19800 

c 

19900 

20000 

c 

20100 

c 

20200 

c 

20300 

20400 

c 

20500 

900 

20600 

901 

20700 

20800 

20900 

21000 

10 

21100 

21200 

21300 

21400 

20 

21500 

21600 

30 

21700 

21800 

21900 

22000 

C 

22100 

C 

22200 

C 

22300 

C 

22400 

C 

22500 

C 

22600 

C 

22700 

22800 

c 

22900 

c 

23000 

c 

23100 

23200 

c 

23300 

c 

23400 

c 

23500 

23600 

23700 

c 

23800 

900 

23900 

24000 

901 

CLOSE CUNIT" 1, Access* 'SEQOUT', FILE* *01, 0*T*) 
CONTINUE 

return 

ENO 


INTERACTIVE  USAGEl 

1)  STORE  Q values  NQ/VES 


NO  additional  routines  used 


***<»stf-**»************* 

subroutine  ODSKROCIDUM) 


READS  U data  from  DISK  FILE  01, DAT 

rouH  * 0 ON  return  for  no  previously  stored  data 
* I ON  return  if  0 values  are  read  from  disk 


common/ A2/0 1 ,02,03, I ST ART, I END, ID J, MO 


USES  COmmqN/AR/  also 


COMMON/ A«/QR (2, 2, 201) ,01(2,2,201) 


F0RMAT(X,II) 

FORMAT(X,Il ,X, It ,X, I3,X,E14.T,X,Et«,7) 

OPEN (UNI T* J , ACCESS* 'SeoiN', file* '01 ,0AT*) 

REAO(l,900)  lOUM 

IF(1OUM.EO,0)  GO  TO  30 

REAU(1,R01,ENO«20)  I,J,N,SR,SI 

QR(I,J,N)*SR 

QI(I,J,N)*S1 

GO  TO  10 

CONTINUE 

IENO*N 

CONTINUE 

Close (UNIT *1, ACCESS* ' sEo IN', file* '01, OAT ») 

return 

End 


NO  interactive  usage 


NO  additional  routines  used 


subroutine  out 


outputs  COMmON/AI/  values 


COMMON/ A2/OmEG0,OMEGNO,OTOmeG, ISTaRT, ieno,np,mo 


USES  COmmon/AI/  also 


C0MM0N/Al/G(2,2,b,2) ,XKB(2,2.3,2) ,XLB(2,2,S,2), 
1 XKA(2,2) ,XLA(2,2) ,F(2,2) 


FORMAT(X,'KA(',H,',',n,')  * ' , FT  . 2 , 3X  , ' L A ( ' , 1 1 , ' , ' , 1 1 , 
')  * ',F7,2,3X,'F(',11,',',II,')  * ',F7,2) 
FORMAT(10X,FT,2,3X,F7,2,3X,F7.2,3x,F7,2,3x,FT,2) 


1 


10 


»< 


:WiV 


K'. » ‘ 


’ / 

iJ 


I 

I 


2aioe 

902 

2«2an 

9t»J 

•J4«*Tttl)t,r7.2,5*.»7,?,3)i,F7,2) 

2«3aa 

9C4 

► :4-*T(K,'«nt',n,',',n.')  • *,32(*-*)) 

24aae 

905 

f:-"»T(x,'U6C*,n.'.',n.')  • *,32('-'>) 

24Saa 

906 

F-.W>1*T  («) 

24(>aa 

9b7 

F’nwMiT  (X,  MIME&P  • ^7.2,  X, '0**EG90  • ' , ^ 7 , 2,  X , » OTOHEC  • • 

247ec> 

/ .P7,2, «, ■ '.Ji.X.'NP  • »,I3) 

24sa? 

9h6 

'Qk“XT  (X, 'The  OHlltB  Of  THE  STSTEH  IS  • '»1D 

24<»'ja 

90ii 

f 09“*T {X , 'Fo«  Output  t»pe:') 

asaap 

9 J l> 

FCB“»T(X,*  0 F04  NONE.l  FOP  TTv.2  F[)»  LP'J 

2siaa 

9J  J 

FOHHiTtx.'TVPET  a TO  C06T1NUE,  TTPCi  1 TO  00  HGXlN'l 

252ae 

912 

f UPHiT  fin 

253a0 

913 

F0W-1XT fx, 'FOR  SySTEh  OUTPUT  TYPE  O.  FOR  0 OUTPUT  TYPE  I'l 

254aa 

5 

*YP£  909 

25Saa 

’yPE  9t0 

2Sbaa 

iCCERT  »,Il 

25700 

fF(It.£O,0)  return 

25800 

»»Pt  906 

25900 

*ypE  913 

26000 

ACCEPT  .,12 

26100 

IF  tla.EO.n  Call  ocutch  J 

26200 

:Fa2.Er..n  go  to  ro 

26300 

n«n-»— 

26400 

iFni.F.Q,n  GO  TO  lao 

26500 

DO  10  1»1,H0 

26600 

00  10  JAl.^O 

26700 

TYPE  906 

26000 

TYPE  90l,6(l,J,5,l},6(I,J,4,n,G(l,J,3,l),G(l,J,2,l),C(l 

26400 

TYPE  902,  I,J 

27000 

TYRE  9ei,G(l,J,5,2),G(I,J,4.2),G(l,J,3,2)>6(I»Ji2,2),G(l 

27100 

TYPE  906 

27200 

10 

CONTINUE 

27300 

00  20  Ip1.no 

27400 

DO  20  Jp1,NQ 

27500 

Type  9o6 

27600 

type  9O0,l,J,IK4(I,2),X,J,XU4niJ>»I»J«F(lf J) 

27700 

20 

CONTINUE 

27000 

DO  30  IPI.NO 

27900 

DO  30  JP1,N0 

20000 

Type  906 

20100 

Type  9e3,IKH(l,J,3,l),IKB(I,J,2,l),X40(t,J,l,l) 

20200 

type  906,1,4 

20300 

type  90S,xk0(1,J,1,2).XKO(X,4,2,2),XKB(I,J,1,2) 

20400 

TYPE  906 

20500 

30 

CONTINUE 

20600 

Do  40  1*1, NO 

20700 

DO  40  4«1,N0 

20000 

Type  906 

20900 

type  903,XlO(I.J,3,1),XLB(1,J,2,1),XLB(I,J,1,1) 

29000 

type  905,1,4 

29100 

Type  903,Xt9(I,4,1.21,XtB(l,4,2,2),XLB(I,4,l,2) 

29200 

type  906 

29300 

40 

CONTINUE 

29400 

45 

Type  907, Onego, onegno.otoneg.ieno.np 

29500 

Type  90b,no 

29600 

90 

continue 

29700 

Type  9u 

29000 

accept  912,I0UN 

29900 

iFdOUH.EO.l)  GO  TO  5 

30000 

IF(I1,EO,0)  0ETU0N 

-iMiL. 
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30100  100 

30200 

30300 

30000 

30S00 

30600 

30700 

30B00 

30900  lie 

31000 

31100 

31200 

31300 

31400  12E 

31500 

31600 

31700 

31600 

31900 

32000 

32100 

32200  ise 

32300 

32400 

32500 

32600 

32700 

32600 

32900 

33000  14e 

33100  14! 

33200 
33300 
33400 
33500 
33600 
33700 
33600  C 
33900  C 
34000  C 
34100  C 
34200  C 
34300  C 
34400  C 
34500  C 
34600  C 
34700  C 
34600 
34900  C 
35000  C 
35100  C 
35200  C 
35300  C 
35400 
35500 

35600  940 

35700 

35800 

35900 

36000 


continue 

OQ  110  I»1|H0 
DO  110 
PRINT  90b 

PRINT  90l(G(IrJr5>l}(G(I«J»4»l)iG(I|J)3(l)rS(IrJ«2it)>G(Iy<7flit) 
PRINT  902,  I,J 

PRINT  901,C(I,J,5,2).G(I,J,4,2),G(I,J.3,2),G(I«J,2,2)*G(I,J,1|2) 

PRINT  906 

CONTINUE 

DO  120  I^liNO 

00  120  J»1»N0 

PRINT  906 

PRINT 

CONTINUE 

00  130  lalfNO 

00  130  Jal,MO 

PRINT  906 

PRINT  903«XKR(I,J,3,1>,XKB(I, J,2,1),XKB(I,J,1>1) 

PRINT  904, I, J 

PRINT  903,XKB(I,J,3,2),XKB(I,J,2|2).XK6(I,J,1,2) 

PRINT  906 
CONTINUE 
DO  140  III, HO 
DO  140  J«1,H0 
PRINT  906 

PRINT  903.XL8(I,J,3,n,XLB(I,J,2,l),XtB(I,J,l>l) 

PRINT  905, 1,J 

PRINT  903,XLB(1,J.3,2),XI.B(I,J,2,2I,XLB(I,J.1,2) 

PRINT  906 
CONTINUE 

PRINT  907,OHEG0,OMEGND,OTOMEG, IEN0,NP 

PRINT  908, MO 

type  911 

4CCEPT  •,I0UH 

lEdOuHTE-Q.l)  GO  TO  5 

return 

END 

Interactive  usage 

1)  OUTPUT  DESIRED  NONE/TTY/LP 

2)  OUTPUT  SYSTEM/0  1 ~ 4»  _ 

3)  00  AGAIN  NO/YES  ^ 

NO  additional  routines  used 


subroutine  finish 


STOPS  EXECUTION  OF  THE  MAIN  PROGRAM 


NO  COMMON  Block  usage 


CALL  ANMOOE 
TYPE  940 

FoRMATCX.'TOU  are  DONEISEST  print  out  SYSTEMIII*) 

CALL  OUT 

CALL  NEWPAC 

CALL  ODSKmR 

CALL  FINITT(0,T0O) 


Sbiei  STOP 

362SS  return 

sssee  eno 

36400  C •*••*«*.*••••••••••••••••••*•****•*•• 

36500  C NO  interactive  USAGE 

36600  C 

36700  C AOOITlONAt  routines  CALLED  I OUT , ODSKMH 

36800  C •.•.,,«,**.,*.*••»*.*•**••***»•*•**•* 
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Table  A. 2. 4.  Autopilot  Design  Program 


00100 

C 

00200 

c 

00300 

c 

00400 

c 

00300 

c 

This  is  The  haih  PR0CR4H  for  the  autopilot  design 

00600 

c 

00T00 

c 

00000 

IMPLICIT  COhPlEX(C) 

00400 

C0MH0N/Al/r.(?,2,6,2)  ,XXB(2,2,3.2)  »XLB(2,2,3,2)  , 

01000 

1 XKA(?,2j,xLA(2,2),F(2,2) 

01100 

COMMON/A2/OHeG0,OMtGND,OTOMEG, ISTARTi IEN0,NP,H0 

01200 

COMMON/ AS j /Cl (2,2) ,C2(2,2) ,31 (2),S2(2) 

01300 

C0MM0N/A32/XW1 (2,2,201) , XU (2,2,201) 

01400 

C0MM0N/ASS/XR2 (2,2,201 ) ,XI2 (2,2,201) 

01500 

C0MH0N/A4/QH (2,2, 201) ,01(2,2,201) 

01600 

C0Hf10N/A5/HAN0(2,2,2ai  ,2) 

01700 

COMHON/A6/nsT(2,20l  ,2) 

01B00 

C0Mm0n/a7/XS(2,2) ,XC(2,2#6,2) 

01400 

c 

02000 

CALL  iMTT(120) 

02100 

CALL  BINTTT 

02200 

Call  0ATA12 

02300 

IENOB(nHEGNO-OMEG0)/OTOHEGtl 

02400 

IF(1ENO.GT,201)  IENOa201 

02500 

CALL  OOSKRO(IOUM) 

02600 

IF  (IDUM.EO.l)  TYPE  BBS 

02T00 

68B 

FORmat(X,'0  values  here  obtained  FROM  DISK') 

02000 

IF(IOUM,eo,l)  60  TO  10 

02400 

DO  5 1*1, HO 

03000 

00  S J«1,H0 

03100 

XS(1,J)«0.0 

03200 

lE(l.EO.J)  XS(I,J)*1,0 

03300 

5 

CONTINUE 

03400 

CALL  HAP(0) 

03500 

CALL  MIMP 

03600 

CALL  OINV 

03700 

10 

CONTINUE 

03000 

TYPE  900 

03900 

900 

F0RMAT(X,'TYPEI  0 TO  STOP') 

04000 

Type  901 

04100 

901 

FORhAT(X,»TYPEI  1 FOR  OPEN  LOOP  PACKAGE*) 

04200 

TYPE  902 

04300 

932 

FO»MAT(X,*TVPEI  2 FOR  INPUT  COMPENSATOR  PACKAGE*) 

04400 

TYPE  903 

04500 

903 

FohhAT(X,*TyPEI  3 FOR  CLOSED  LOOP  PACKAGE*) 

04600 

type  904 

04700 

904 

FoRmAT(X, *TyPEI  4 FOR  OUTPUT*) 

04000 

TYPE  905 

04900 

905 

FORmAT(X, 'TyPEI  5 TO  modify  INTEGER  C0MM0N/A2/  VALUES*) 

05000 

ACCEPT  »,inuH 

05100 

IF(IOuH.EO,0)  CALL  FINISH 

05200 

IFdnuM.EQ.l)  CALL  OLISP 

05300 

1F(1DUM,E0.2)  LALL  IFOCP 

05400 

IF(IOu«,£Q,S)  call  CLBP 

05500 

iF(inuH,eo.4)  CALL  OUT 

05600 

1F(10UH.E0.5)  CALL  M001A2 

05700 

60  TO  10 

0S000 

STOP 

05900 

END 

00000 

C 

A 


0010a 

00200 

C 

5U8ff0u7INE  OLISP 

OOS00 

00400 

c 

1MPLICI7  COHPtEXCC) 

O0SOO 

C0HHnN/Al/r.(2,2,6,2) .XKU(2.2.3,2),XUB(2|2|S,2)» 

00400 

1 XK4{2,2).XLA(2,2),7(2,2) 

00700 

COHMON/*2/OM£G0,OHtU6O,OTOHEC,lST*0T»I6NO,NP,HO 

00800 

C0MM0n/*31/C1(2,?),C2(2,2),S1(2),S2(2) 

00400 

common/ A32/XH1 (2.2.201), XII (2,2,201) 

01000 

CUMhOn/aS3/XR2(2,2,201),XI2(2,2,2O1) 

01100 

COMMON/ A4/QP (2, 2, 201), 01 (2, 2, 201) 

01200 

CQMMaN/45/HAN0(2,2,2ei,2) 

01300 

COMMON/ *6/057(2,20 1,2) 

01400 

O1S00 

c 

CUMM0n/*7/XS(2,2),XC(2,2,6,2) 

01600 

DO  110  I«1,M0 

01700 

DO  110  Jal.MO 

01800 

XS(I,J)*XK*(X,J) 

01400 

DO  100  N»1,1EN0 

02000 

*«l(l,J,N)»OP(I,J,N) 

02100 

Xll(I,J,N)»OI(X,J,N) 

02200 

100 

continue 

02300 

110 

C0N7INUE 

02400 

C*l.t  CERSN 

02500 

type  410 

02600 

410 

Format (X, 'DEALING  W/  OPEN  LOOP  INPUT  CONST,  COHP*) 

02700 

120 

CONTINUE 

02800 

type  400 

02400 

400 

FORMATIX.^TYPEI  0 TO  RETURN,  1 TO  PLOT,  2 TO  HOOXFY*) 

03000 

ACCEPT  •,10UH 

03100 

CALL  NEwPAG 

03200 

XFIXDUM.EQ.O)  GO  TO  140 

03300 

IF(IOUM.EO.l)  GO  TO  130 

03400 

CALL  HOOXS 

03500 

Go  TO  120 

03600 

130 

CONTINUE 

03700 

CALL  HXMP 

03800 

CALL  GERSH 

03400 

CALL  ORNI 

04000 

GO  TO  120 

04100 

140 

CONTINUE 

04200 

00  150  1>1,H0 

04300 

00  150  J*1,M0 

04400 

XK*(X,J)>XS(X,J) 

04500 

150 

Continue 

04600 

144 

CONTINUE 

04700 

return 

04S00 

END 
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00100 

c 

00200 

subroutine  cusp 

00300 

c 

00100 

IHPlICIT  COmPuEx(C) 

00500 

COHhon/41/G(2,2,6,2) »XKB(2«2.3,2) ,XU0(2«2,3,2)i 

00600 

1 XK4(2,2),xL*(?r2),F(2,2) 

00700 

COHNON/A2/OMeG0,nMEGNO,t)TOMES,IST*RT, lENOiNPjMO 

00000 

C0MM0N/*3|/Ct (2,?) ,C2C2,2) tSl(2] |S2(2) 

00900 

C0NM0N/A32/XR1 (a,2,201),XJl (2,2,201) 

01000 

CONMON/A33/XR2(2,2,201),XI2(2,2,201) 

01100 

common/ A4/0H (2, 2, 201), 01 (2, 2, 201) 

01200 

COMHON/*5/8and(2,2,201,2) 

01300 

common/ A6/0ST (2,201,2) 

01400 

C0MM0n/A7/XS(2,2),XC(2,2,6,2) 

01500 

c 

01600 

00  310  I«1,M0 

01700 

DO  310  J*1,H0 

01000 

XS(1,J)«P(I,J) 

01900 

310 

CONTINUE 

02000 

CALL  GERSN 

02100 

CAUL  05TR0W 

02200 

TYPE  930 

02300 

930 

F0RMAT(X, 'DEALING  M/  CLOSED  LOOP  SYSTEM*) 

02400 

320 

CONTINUE 

02500 

type  900 

02600 

900 

PORMAT(X,'TYPEI  0 TO  RETURN,  1 TO  PLOT,  2 TO  MOOIPT*) 

02700 

ACCEPT  *,  lOUM 

02000 

CALL  NEwPAG 

02900 

lEdDUM.EQ.O)  60  TO  340 

03000 

IFdOUM.EO.l)  GO  TO  330 

03100 

CALL  MODXS 

03200 

Go  TO  320 

03300 

330 

CONTINUE 

03400 

CALL  OSTRQM 

03500 

CALL  PTSIGl 

03600 

60  TO  320 

03700 

340 

CONTINUE 

03000 

DO  350  l>t,MO 

03900 

00  350  J*1,M0 

04000 

F(I,J)-XS(I,J) 

04100 

350 

CONTINUE 

04200 

399 

Continue 

04300 

RETURN 

04400 

END 

110 


eeiee 

B020e 

C 

subroutine  0*7*12 

eeaee 

w 

C0HN0N/*1/G(2,2|6,2) •XKB(2»2.3,2),X|.6(2,2,3,2}i 

00S0B 

1 *K*(2,2),*L*t2,2),F(2,2) 

00b00 

common/ *2/OnEG0, OHEGNO,OTOHEG, IST*RT,IENO, NP, HO 

00700 

c 

SETS  COMMON  BLOCK  /*2/  V*LUES 

00600 

OMEG0S0.0  . 

00900 

UMEGN0*0,16 

01000 

DTUMEG«0,0006 

01100 

IST*RT»1 

01200 

IENO«201 

01300 

NP»10 

01400 

H0a2 

01S00 

c 

SETS  *LL  DEMON,  OF  G(S)«CH*R,  POLYNOHI*L 

01600 

00  10  lalfHO 

01700 

DO  10  3*1, HO 

01600 

G(1»J,1,2)«0.0 

01900 

G(I, J,2,2]«-.00166 

02000 

6(I,J,3,2)«,365 

02100 

6(I,J,*,2)*.364 

02200 

G(I,J,5,2).,829 

02300 

G(I,J,6,2)»1.0 

02400 

10 

continue 

02500 

c 

SETS  1N0IVIDU*L  NUMER*T0R  ENTRIES  OF  G(S) 

02600 

G(l,l,l,n«-,0195 

02700 

cn,i,2,n«-.00893 

02600 

G(l,i,3,n»-,3121 

02900 

C(l,l,4,l)«-,3807 

03000 

e£l,l,5,l)»0.0 

03100 

GC1»1,6,1)»0.0 

03200 

G(l,2,l,l)«0.0 

03300 

G(1,2,2,1)«-.4*7 

03400 

S(I,2,3,1)»-.0503 

03500 

Ca,2,4,l)»-.0404 

03600 

6(l,2,5,na0,0 

03700 

G(1,2,6,I)»0,0 

03600 

G(2,l,l,l)». 022635 

03900 

Gt2, 1,2,1)*-. 000566 

04000 

G(2,],3,na-.017S 

04100 

C{2,l,4,n».067l5 

04200 

G(2,l,5,l)»0,0 

04300 

G(2,i,6,1)«0,0 

04400 

G(2,2,l,l)«0,0 

04500 

O{2,2,2,l)«.b02 

04600 

G(2,2,3,t)a.06237 

04700 

G(2,2,4,nai,567 

04600 

C(2,2,5,l)»0,0 

04900 

G(2,2,6,l}a0.0 

05000 

c 

SETS  XK*at 

05100 

KK*(1, 1)«1.0 

05200 

XK*(1,2)«0.0 

05300 

lK*(2,l)ae,o 

05400 

VK*(2,2)*1.0 

05500 

RETURN 

05600 

ENO 

05700 

c 

I 


00100  c 
00200 
00300  C 
00000  C 
00500  C 
00600  C 
00T00 
00000 
00000 
01000 
01100  C 
01200 
01300 
01400 
01500 

01600  10 
01700 
01000 
01400 
02000 
02100 
02200 

02300  900 

02000 
02500 
02600 
02700 
02600 

02900  IS 

03000 

03100 

03200 

03300 

03000  20 

03500  30 

03600 

03700 

03600  C 


subroutine  OINV 


TAKES  'LONG  HANO'  INVERSE  OF  0 
FOR  A 2 * 2 SYSTEM  ONLY 


implicit  COmPlEX(C) 

COMMON/A2/OmEG0,OmEGNO«OTOMEG,ISTART* IENO,NP,MO 
COMMON/ a3i /Cl (2,2) ,C2(2,2) ,si(2) ,S2(2) 

COMMON/ A0/OW (2, 2, 201) ,01(2,2,201) 

00  30  Nal,lENO 

00  10  Ial,2 
DO  10  Jil,2 

Cl (I,J)«CMPLV(QR(1,J,N),01(X,J,N)) 

CONTINUE 

XREaL«0R(1,1,N)*QR(2,2,N)4'01(1,2,N) aOI(2,1,N) 

1 -Uld ,1,N)*0I (2,2,N).QR(1,2,N)*QR(2,1,N) 
XIMAG*0I(l,l,N)*0R(2,2,N)A0R(l,l,N)*01(2,a,N) 

1 •0R(l,2,N)*0I(2,l,N).QI(l,2,N)*0R(2,i,N) 
CCC«CMPLX(XNEaL,X1MAG) 

IF  (CCC.EO.CMPL*(0. 0,0.0))  TYPE  900, N 

F0RMAT(X,'Q  singular  upon  INVEPSION  P ',I3) 

IF(CCC.EO.CmPlX(0, 0,0,0))  GO  TO  IS 

C2(l,l)«Cl (2,2)/CCC 

C2(1,2)*-C1(1,2)/CCC 

C2(2, 1)»-C1 (2,1)/CCC 

C2(2,2)«Cl(l,l)/CCC 

CONTINUE 

00  20  I«l,2 

00  20  J«l,2 

0P(I,J,N)aREAL(C2(X,J)) 

0I(I,J,N)«AIMAC(C2(1,J)) 

continue 

CONTINUE 

RETUHN 

end 


